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Abstract 

Monotone lattice recurrence relations such as the Frenkel-Kontorova lattice, arise in 
Hamiltonian lattice mechanics, as models for ferromagnetism and as discretization of elliptic 
PDEs. Mathematically, they are a multi-dimensional counterpart of monotone twist maps. 

Such recurrence relations often admit a variational structure, so that the solutions x : 
Z'* — >■ R are the stationary points of a formal action function W{x). Given any rotation 
vector tj G R'', classical Aubry-Mather theory establishes the existence of a large collection 
of solutions of VW(x) = of rotation vector cj. For irrational lo, this is the well-known 
Aubry-Mather set. It consists of global minimizers and it may have gaps. 

In this paper, we study the parabolic gradient flow ^ = — VVK(a;) and we will prove that 
every Aubry-Mather set can be interpolated by a continuous gradient-flow invariant family, 
the so-called 'ghost circle'. The existence of these ghost circles is known in dimension d = 1, 
for rational rotation vectors and Morse action functions. The main technical result of this 
paper is therefore a compactness theorem for lattice ghost circles, based on a parabolic 
Harnack inequality for the gradient flow. This implies the existence of lattice ghost circles 
of arbitrary rotation vectors and for arbitrary actions. 

As a consequence, we can give a simple proof of the fact that when an Aubry-Mather 
set has a gap, then this gap must be filled with minimizers, or contain a non-minimizing 
solution. 



1 Introduction and outline 

In this paper we are interested in variational monotone lattice recurrence relations. Before 
introducing such recurrence relations in full generality, let us discuss as an example the 
so-called d-dimensional Frenkel-Kontorova lattice. Here, the goal is to find a d-dimensional 
"lattice configuration" a; : Z'' — >■ R that satisfies 

V'{x,) - (Ai), = for all i e Z"*. (1.1) 

In the equation above, the smooth function : R — >■ R satisfies + 1) = V{^) for all 
^ £ R. It has the interpretation of a periodic onsite potential. Setting ||i|| := X]fc=i ^he 
discrete Laplace operator A : R^ — >■ R^ is defined as 

(Ax), := ^ [xj ~ x^) for aU i £ Z'*. (1.2) 

l|j-i||=i 

One could think of equation as a naive discretization of the nonlinear elliptic partial 

differential equation V' {u) — Au = for a function m : R'' — >■ R and Xi = u{i). 

At the same time, equation is relevant for statistical mechanics, because it is related 
to the Frenkel-Kontorova Hamiltonian lattice differential equation 

^ + V'{x,) - (Ax), = for all i € Z". (1.3) 
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This differential equation describes the motion of particles under the competing influence of 
an onsite periodic potential field and nearest neighbor attraction. Obviously, equation 
describes its stationary solutions. 

Finally, in dimension d — 1, the solutions of equation ( |1.1[ ) correspond to orbits of the 
famous Chirikov standard map Tv of the annulus. This correspondence is explained in some 
detail in Appendix A. 

The Frenkel-Kontorova problem is an example from a quite general class of lattice 

recurrence relations to which the results of this paper apply. These are recurrence relations 
for which there exists, for every j G Z'', a real-valued "local potential" function Sj : — > R 
so that the relation can be written in the form 

^ ^^^(x) = for alH G Z''. (1.4) 

It turns out that for the Frenkel-Kontorova problem such local potentials exist and it 

is easy to check that they are given by 



5*. 



{x):=V{x,) + l^ J2 i^k-x,f. (1.5) 



l|fc-j||=i 



For the general problem ( |1.4[ ), the functions Sj{x) will be required to satisfy some rather 
restrictive hypotheses that will be explained in detail in Section [2] Physically, the most 
important of these hypotheses is the monotonicity condition. It is a discrete analogue of 
ellipticity for a PDE. Among the more technical hypotheses is one that guarantees that the 



sums in expression ( 1.4| are finite. For the purpose of this introduction, it probably suffices 



to say that the potentials ( 1.5 1 of Frenkel-Kontorova are prototypical for the Sj{x) that we 
have in mind. 



It is important to observe that the solutions of (1.4 1 are precisely the stationary points 
of the formal sum 

W{x) := ^ SAx). (1.6) 
This follows because differentiation of (|1.6[ ) with respect to Xi produces exactly equation 



(1.41 and it explains why solutions to (1.4 1 are sometimes called stationary configurations. 



In the case that the periodic onsite potential V^(^) vanishes, the Frenkel-Kontorova equation 
( |1.1[ ) reduces to the discrete Laplace equation Ax = 0, for which it is easy to point out 
solutions. For instance, when ^ G R is an arbitrary number and cj G R'' is an arbitrary 
vector, then the linear functions x^'^ : Z"* — >■ R defined by 

:=? + (c.,i) 

obviously satisfy Ax = 0. It moreover turns out that the x'^'^ are action-minimizers, in the 
sense that for every finite subset B C Z'' and every y : Z"* — >■ R with support in B, it holds 
that 

^ (s,(x-'« + y)-5,(:r-'«)) >0. 

Note that this sum is actually finite and can be interpreted as W{x^'^ + y) ~ W(x^'^). 

Definition 1.1. Let a:: : Z'' — >■ R be a d-dimensional configuration. We say that tj G R'' is 
the rotation vector of x if for all i , the limit 

lim exists and is equal to {ijJ,i). 
Clearly, the rotation vector of j:"'* is equal to lo. On the other hand, in dimension d 7^ 1, a 



solution to ( 1.1 1 does not necessarily have a rotation vector. An example is the hyperbolic 
configuration x defined by x'l = i\i2 ■ ■ ■ id-iid which solves Ax = 0. 

In Aubry-Mather theory, one is interested, among others, in answering the following ques- 
tions: given a collection of local potentials Sj (x) satisfying the assumptions of Section |2j a 
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number ^ G R and a vector G R , does there always exist a solution x to equation (1.4 1 
with rotation vector u and initial condition xo = ^? And if so, what is the structure of the 
solution set? 

A rather complete answer to these questions is known. It turns out that solutions to 
11.4 1 of all rotation vectors cj £ R'' exist. For example, it was shown by Bangert that 
when a; e R''\Q'' is irrational, then there exists a unique nonempty collection of "recurrent" 
action-minimizers of rotation vector uj. This is the Aubry-Mather set of rotation vector uj. 
It is totally ordered, but may contain "gaps". That is, given an arbitrary ^ G R, it may 
happen that the Aubry-Mather set of rotation vector uj does not contain any configuration 
X satisfying the initial condition xq = It is known that in this case, the Aubry-Mather 
set is actually a Cantor set. 

The basics of this classical theory will be reviewed in Sections |3] and |4] of this paper. In 
Sectionjsj we will study Birkhoff configurations, examples of which are the action-minimizing 
configurations of the Aubry-Mather sets. In Section [3. 2| we will moreover prove some new 
results for d-dimensional periodic Birkhoff configurations. In Section |4j we will examine 
minimizing configurations and for completeness, we will reprove the classical result that 
global minimizers of every rotation vector exist and we will examine the properties of the 
Aubry-Mather set. 

To investigate the existence of stationary configurations in the gaps of the Aubry-Mather 
sets, we propose to study the gradient fiow of the formal action function, i.e. the fiow of the 
differential equation 

^ = ~VW{x) . 
dt ^ ' 

It was shown by Gole [7] that this fiow is well-defined on a suitable subspace X C R^ of 
configurations that contains all Birkhoff configurations. We will prove some regularity results 
for the gradient flow in Section [S] and we will discuss some of its qualitative properties in 
Section [6] The most notable of these is a strong monotonicity property or strong parabolic 
comparison principle, see Theorem |6.2[ 

The principal goal of this paper is then to prove the existence of a continuous one- 
dimensional gradient-fiow invariant family of configurations that contains the Aubry-Mather 
set of rotation vector uj. Such an interpolating family will be called a ghost circle and denoted 

C R^ . The precise definition of a ghost circle is given in Section [t] 

Ghost circles were already constructed for twist maps by Gole [8]. Hence, they are 
well-known to exist in dimension d = 1. Gole starts his construction by assuming that 
= I G Q is rational and that an appropriate periodic action function W-p^q{x) is a Morse 
function. Under these assumptions, the existence of a periodic ghost circle follows from a 
combination of topological arguments and the parabolic comparison principle of the gradient 
flow. In Section [8.2[ we will imitate the construction of these periodic Morse ghost circles in 
dimension d 7^ 1. Each of these periodic ghost circles contains at least one global minimizer. 

Our first main result is contained in Section |8.1| It generalizes results of Gole on 
twist maps and it roughly states that for every rational bj G Q'' and for every collection 
of potentials Sj{x), one can find arbitrarily small perturbations of the Sj{x) that turn the 
periodic action Wp^q{x) into a Morse function. This statement is nontrivial in dimension 
d 7^ 1 and it holds because of group theoretic reasons that will explained in Section |3.2[ 

The most important technical result of this paper is nevertheless a compactness theorem 
for ghost circles. It is presented in Section [9] It says that when the rotation vectors u)„ 
converge to a rotation vector ujoa and the local potentials S" converge to potentials Sj° and 
there exist ghost circles r„ for the potentials of rotation vector aj„, then there is a ghost 
circle Too for the potentials Sj° of rotation vector ljoo ■ Moreover, a subsequence of the r„ 
actually converges to Too in a sense to be made precise. Together, all of the above shows 
that there are ghost circles of every rotation vector and for arbitrary potentials. Again, they 
contain at least one minimizer and hence the entire Aubry-Mather set of rotation vector oj. 

A similar compactness result was proved by Gole, see 0, for twist maps. The proof of 
this "monotone convergence theorem for ghost circles" relies on the fact that over time, two 
different solutions of the gradient flow must decrease their number of intersections. Hence, 
this proof is purely one-dimensional. Our proof, on the other hand, only depends on a 
quantitative version of the parabolic comparison principle, a so-called Harnack inequality. 
This inequality is stated and proved in Theorem |6.4| 

As a consequence, we show in Section \To\ that when the Aubry-Mather set is a Cantor 
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set, then its gaps must either be completely foliated by minimizers, or contain at least one 
non-minimizing solution to ( |1.4[ ). 
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2 Problem setup 

Let us at this point introduce the generalized Frenkel-Kontorova lattice recurrence relations 
that we want to consider in this paper. 

As was discussed before, we will assume that for all j £ Z"* there is a function Sj that 
assigns a real value to every d-dimensional configuration: 



These functions are required to have the conditions A-E described below and are called local 
potentials. 

To formulate the first condition, let us assume that a finite subset _B C Z'* and an m > 
times continuously differentiable function s : — >■ R are given. Then we can define a 
function S : R^ — >■ R by setting S{x) :— s{x\b)- This is just a way of saying that S depends 
only on the finitely many variables Xi for which i £ B. 

Such an S has some convenient properties, most notably that it is continuous in the 
topology of pointwise convergence: if x",x°° G R^ and lim„^cx) 2:" ~ x°° pointwise, then 
obviously also lim„^oo 5'(a;") — S{x°°). 

Moreover, it makes sense to speak of the partial derivatives of the function S: if ji , . . . , jk G 
Z'*, with < < m, is a collection of lattice points, then the partial derivative djj^,....j^S : 
R^ — >■ R can simply be defined as 



dji,...,jkS{x) 



(9ji,...,jfcS)(a^|s) if jl, ■ ■ -Jk e B, 
otherwise. 



These partial derivatives are also continuous with respect to pointwise convergence 

lt=i I**:!' f*-"^ i G Z'', and define Bj 



Finally, we recall the definition ||i|| := X]fc=i Nfcli ^'^^ * ^ ^""i ^^'^ define Bj := {k G 



Z^ 1 \\k^j\\<r}. 

With all this in mind, we can formulate our first condition. 

A. The functions Sj are twice continuously differentiable and of finite range. That is, 
there is an < r < 00 and for every j G Z"* there is a twice continuously differentiable 
function Sj : R"^j — >■ R such that Sj{x) = Sj{x\B':)- 

In other words, the function Sj depends only on the finitely many variables Xk with | j 1 1 < 
r. Hence, Sj (x) has the interpretation of the "local energy" of the configuration x at lattice 
site j and we think of r as the finite range of the interaction. 

To formulate condition B, it is convenient to introduce an action of Z'' x Z on R^ by 
"shifts": 

Definition 2.1. Let G Z'* and / G Z. Then the shift operator Tk,i ■ R^'' -i^ R^*" is defined 
by 

{Tk,ix)i ~ Xi+k + I ■ 

Clearly, the graph of Tk,ix, viewed as a subset of Z** x R, is obtained by shifting the graph 
of X over the integer vector {—k,l). This explains why the Tk^i are called shift operators. 

B. The functions Sj are shift-invariant: Sj{Tk.ix) = Sj+k{x) for all j, k and I. 

In fact, invariance of the Sj under ro,i just means that Sj{x) — Sj{x + l^d) for all j, which 
means that Sj descends to a function on R^ /Z. Invariance of the Sj under the shifts Tkfi 
expresses the maximal spatial homogeneity of the local potentials. In fact, once one of the 
Sj is given, for instance So, then all the others are determined. 
The next condition ensures the growth of the 5*^ at infinity: 
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C. The functions Sj are bounded from below and coercive in the following sense: for all 
k with I |fc — J 1 1 — 1, 

lim Sj (x) — oo . 

\Xf^—Xj I — >oo 

Condition C says that every function x i— >■ Sj{x) is as coercive as it can possibly be under 
the restriction that it satisfies the periodicity condition Sj(ro,ia;) = Sj{x). 
The following condition D is the most essential one: 

D. The functions Sj satisfy the so-called monotonicity condition: 

di^kSj < for all j and all i 7^ fc, while di^kSi < for all | |i — fc| | = 1 . 

Condition D is also called a twist condition or ferromagnetic condition. It says that all mixed 
derivatives of the local potentials are non-positive, while some of them are strictly negative. 
For technical reasons we will also assume: 

E. The Sj have uniformly bounded second derivatives: there is a constant C such that 

\di,kSj \ < C for all i.,j,k. 

As in Section [1] we can now look for stationary configurations corresponding to these po- 
tentials. 



Definition 2.2. A configuration called a stationary point for the local 

potentials Sj if for every finite subset B G Z'^ and every configuration y with support in its 
r-interior := {i e B \ B[ C B}, it holds that 

WB{x + ey) = 0, 



d_ 
de 



where Ws : K R is defined as 

VKs(x-):=^5,(x). (2.7) 

In fact, by difi^erentiating Wb with respect to an Xi with i G B^'^\ one obtains that x is 
a stationary point for the Sj if and only if it satisfies the variational monotone recurrence 
relation 

^ ftS'j (x) = for alH e Z''. (2.8) 

\\j — i\\<r 



The goal of this paper is to find solutions of (2.81 and while doing so, we will exploit the 
variational principle that underlies it. 

By the way, (12. Si is called monotone because condition C guarantees that the derivative 



of the left hand side of (2.81 with respect to any of the Xk with fc 7^ i, is non-positive, while 
it is strictly negative if ||fc — ijj = 1. 

Definition |2.2| moreover inspires the definition of a special type of solutions to ( |2.8[ ): 

Definition 2.3. A configuration a; : Z'' — >■ R is called a global minimizer or ground state for 
the potentials Sj if for every finite subset B C Z'' and every j/ : Z'' — >■ R with support in 

Wb{x + y) - Wb{x) = (^^(^ + V)' 5^ ^ ■ 

Clearly, global minimizers are automatically stationary and hence satisfy the recurrence 
relation ( 2.8 1. 



Example 2.4. It is easy to check that the Frenkel-Kontorova potentials given in ( 1.5 1 satisfy 
conditions A-E. In fact, the range of interaction is r = 1, and dj^kSj — — for \\j — k\\ = 1. 

In the particular case that V{^) = all solutions of ( |1.1| are actually global minimizers. 
This follows because every y 1— > Wb{x + y) is strictly convex if V{£) = and hence only has 
one stationary point, which is minimizing. 
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3 Spaces of configurations 



In this section, we introduce certain spaces of configurations that are often encountered in 
classical Aubry-Mather theory. We will moreover study some of their properties. Most of 
the definitions and results in this section are standard, but to the best of our knowledge 
Lemma |3.10| and Theorem |3.12| in Section [3.2| are new. We start by recalling the following 
definition: 

Definition 3.1. Let i : Z** — >■ R be a d-dimensional configuration. We say that to — u]{x) £ 
R"* is the rotation vector of x if for all i £ Z'', the limit 

lim exists and is equal to {uj,i). 

n— >oo 72 

The space of configurations with rotation vector u) is denoted 

X„ — {x : Z'' R I uj{x) = ui} . 

3.1 Birkhoff configurations 

We will now introduce the concept of a well-ordered lattice configuration. 

Definition 3.2. On the configuration space R^ we define the relations <, < and <^ by 

• 2; < J/ if a:;^ < for every i € Z'^ . 

• X < y if X < y, but x y. 

• X <^ y li Xi < yi for every i £ U^. 
Similarly for >, > and 

Recall the definition of the shift operators Tk,i : R^ — )■ R^ . The partial orderings defined 
above, now allow us to make the following definition, as in for instance [3] and [11) . 

Definition 3.3. A configuration a; £ R^ is called a Birkhoff configuration or a well-ordered 
configuration, if for all A: £ Z'' and I £ Z, 

either Tk.ix > x or Tk.ix < x. (3.9) 

Definition |3.3| says that the graph of a Birkhoff configuration x does not cross any of its 
integer translates. The space of Birkhoff configurations will be denoted ^? C R^ and it 
inherits the topology of pointwise convergence. Birkhoff configurations will play an essential 
role in the remainder of this paper. Birkhoff configurations of every rotation vector exist: 
for every tj £ R'' the linear configuration x^ defined by xf :— {uj,i) is an example. 

Remarlc 3.4. When : R/Z — >■ R/Z is an orientation preserving circle homeomorphism, 
then it admits a lift to a strictly increasing map iJ" : R — >■ R that satisfies H{(, + 1) — H{^) + 1 
and i7(0modl = /i(^modl). 

Let us now denote by x{^) : Z R the Jf-orbit of ^ £ R, defined by x{^), ■- H'{^). 
Then it is clear that for ^1,^2 £ R with ^1 < ^2, one has that x(^i) <^ x{£^2)- In turn this 
implies that each x{^) is a Birkhoff sequence. Thus, ordering is a very natural concept in 
the theory of circle homeomorphisms. 

The following result is folklore and it goes back to Poincare, who proved it in the case d = 1 
and in the context of circle homeomorphisms, for which it implies that circle homeomor- 
phisms have a unique rotation number. 

For d = 1, the proof of Lemma [3.5| can be found for instance in [S]. For completeness, we 
include the proof for d > 1 here. Lemma [3. 5| says that the graph of a Birkhoff configuration 
X lies uniformly close to the graph of the affine configuration i t-^ xq + {oj , i) ■ 

Lemma 3.5. LetxGM.^ be a Birkhoff configuration. Then x has a rotation vector uj = lo(x) 
and 

\xi - xq- {uj{x),i)\ < 1 . 

Moreover, the map x 1—^ ^{x), B M.'' is continuous with respect to pointwise convergence. 
We write 

B^~{x(£B\ uj{x) = uj} . 
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Proof. We will assume that the result is true for d = 1 and we choose i,j £ Z . Then the 
sequence n i— >■ x„i+j is a one-dimensional BirkhofF sequence and hence its rotation number 
LUij exists and is equal to lim„^oo . Moreover, \xni+j —Xj — {uJij ,n)\ < 1. We first of all 

remark that ujij does not depend on j, and hence can be denoted uji. This follows because 
the Birkhoff property of x ensures that the sequences n i— >■ Xni+j = {Tjfix)ni and n i— > Xni 
do not cross. Now denote by ei = (1, 0, . . . , 0), 62 = (0, 1, 0, . . . , 0), etc. the standard basis 
of and define uj := (ojei, . . • ji^e^)- Then, 

\xi — xo — = k{n,i2,-.-,id) ~ a;(o,i2,...,i^)— iiojei + . . . + a;(o,...,o,ij) — xo — idi^e^\ < 

|a;(ii,i2,...,id) ~ a;(o,i2,...,i^) — iiu)ei\ + . . . + |a::(o,...,o,jd) — ^0 — idi^ej_\ < d . 

This clearly implies that lim„_>oo = while the Birkhoff property of the sequence 

n I— >■ x„i then implies that in fact, \xi — xo — {uj, i)| < 1. 

The continuity of a; i-> a; (a;) follows immediately from the continuity in the one-dimensional 
case. □ 

The following proposition is equally standard. In particular, it will allow us to take limits 
of Birkhoff configurations with rational rotation vectors in order to produce Birkhoff config- 
urations with irrational rotation vectors. 

Recall the action ro_i : s x + 1 on . It can be used to identify sequences that differ 
by an integer. The quotient space is denoted /Z. Note that every element [x] in this 
quotient space has a unique representative x with xo G [0, 1). 

Proposition 3.6. Let K C M.'' be compact and let Bk '■= Uc^gA' Then Bk /Z is compact 
in the topology of pointwise convergence. 

Proof. By definition, B is closed in the topology of pointwise convergence. Moreover, by 



Proposition 3.5 Bk/Z is a closed subset of 



{[x] e R^ /Z \ Xk = xo + k ■ oj + yk with {[x]o,uj,yk) G R/Z x A' x [-1,1]^ } , 

which is compact in the topology of pointwise convergence. This follows from Tychonov's 
theorem. □ 

The following corollary of the compactness of Bk/Z is trivial, but it has important implica- 
tions. 

Corollary 3.7. Let A C R'' fee compact and let B C Z"^ be a finite subset. Assume that 
s : M.^ /Z — > R is a continuous function. Then S : R^ /Z — >■ R defined by S{x) — s(x\b) 
attains its maximum and minimum values on Bk/Z. 

Proof. This follows because such a 5 is continuous with respect to pointwise convergence 
and Bk/Z is compact. □ 

Applied to 5* = di^kSi with ||z — fc|| — 1, and recalling the twist condition di^kSi < 0, 
Corollary |3.7| implies that there is a A > such that di^kSiix) < —A < for all x £ Bk- In 
other words, the twist condition D is automatically uniform on Bk. Similarly, even if one 
does not impose condition E, there is a constant C > such that \di,kSj\ < C for all i,k 
and j, uniformly on Bk- 

We finish this section with a simple and well-known proposition that expresses that the 
number theoretical properties of the rotation vector a; of a Birkhoff configuration x decide 
to a large extent wether Tk^ix > x or Tk,ix < x. 

Proposition 3.8. Let tj £ R'* and x € B^. If {i-o,k) + I > 0, then Tk.ix > x and if 
{uj, k) + I < 0, then Tk.ix < x. 

Proof. Denote by x'^ £ Bu, the linear configuration defined by xf :— Then Tkjx'^ — 

a;" = {k,uj) + I. Suppose for instance that {k,uj) + I > 0, that is that rk,ix'^ ^ a;", but 
assume on the other hand that Tkjx < x. This means that Tk^iX — x < and hence also 
that ri iX — Tk,ix = Tk,o{rk.ix — x) <Q. Thus, rliX — x = {rliX — Tk,ix) + (rk,ix — x) < 0, 
i.e. tI ix < X. By induction we then find that ja; < x, for every n > 1. On the other 
hand, Tkix'^ = x" + n{{k,uj) + I). This contradicts the fact that sup, |r^;(a;" — x)i\ = 
supj |(a;" — x)i\ < \xo\ + 2 is uniformly bounded in n. □ 
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3.2 Periodicity 



It turns out convenient to consider periodic configurations. To define these, let pi, . . . ,Pc £ 
Z'* be < c < d linearly independent integer vectors and let qi,. . . ,qc € Z be c integers. 
Then we set 

Xp,, := {x : Z"* R I Tp^^q^x = a; for all j = 1, . . . , c } . 

We say that a configuration periodic with periods (pi, gi ),..., (pc, ?c). The 

collection of periods of Xp,, is a lattice of rank c, that we denote by 



'/ X 



An element of Xp,, can have a rotation vector, but this rotation vector can not be arbitrary: 
when a; : Z'' — >^ R is a configuration of rotation vector oj and Tpj x = x, then Xnpj = xo—nqj , 
so that lim„_Kx) -^p- = —qj, that is 

{ui,pj) +qj =0 when Xp,q n X;^ . 

Another way to express this is that when Xp,g PI X^^ ^ 0, then Jp,g C /w, where the lattice 
lui is defined as 

lu ■= {{k,l) e Z'^ X Z I {k,uj)+l = 0} C Z'^ X z . 

On the other hand, when x has rotation vector w and (w, k}+l = 0, then this does not imply 
that Tk,ix = x. We therefore define 

X„ := {x £ X„ I Tk^ix — x when (a;, fc) + Z = } . 

The elements of X^ are called maximally periodic as they have all the periods that an element 
of "Ku, can possibly have. X„ is nonempty because it contains the linear configuration a" 
defined by xf = {uj, i). This is true because Tk,ix" = x'^ + (u), k) + I. 

Definition 3.9. A Z-basis (pi, gi), . . . , (pc, 5c) of I,^ is called a collection of principal periods 
for uj. That is, (pi,gi), . . . , (pc,gc) are principal periods for w if and only if Xp,q flXo; = X^,,. 

Of course, a set of principal periods for w G M'* always exists, but it is not unique. 

At this point, let us make some group theoretic remarks. First of all, we remind the 
reader that we can think of the shift operators Tk,i as defining a group action of Z*^ x Z on 
the space of configurations: 

T : (Z'' X Z) X R^"* ^R^**, {{k,l),x) ^ Tk,ix . 

Clearly, because Z'' x Z is Abelian, when Tp .q x = x, then also Tp, ^q.ijk.ix) = t^^ix, and 
thus T leaves Xp,q invariant. Moreover, because the elements of Jp^q fix all elements of Xp,q, 
we have that when x € Xp^q and (fc, I) — {K, L) + J]j "^ifei Ij) for certain integers mj, then 
Tk,ix — tk,lx. This shows that t induces an action of (Z*^ x 1)/ Jp^q on Xp,,. We recall that 
this action is called free if for every (fc, I) ^ Jp,q and every x £ Xp,q it holds that Tk,ix ^ x. 
We now have the following quite obvious characterization of X^, : 

Lemma 3.10. Assume that Xp,, nX;^ ^ 0. Then the T-action of (Z'* x Z)/Jp,q on Xp^q (IX^ 
is free if and only if the {pj,qj) are principal periods for lv , i.e. if and only ifXp^qilX^i — X„. 

Proof. Let us start by assuming that the {pj,qj) are principal periods for w, that is that 
Jp,q = luj- We want to show that then the action of (Z'* x Z)//;^ on Xp,g fl X„ is free. 
But a nontrivial equivalence class in {U^ x Z)//^^ is represented by an element {k,l) with 
{ijj,k) -\- I =^ Q and it is clear that this inequality implies that Tk,ix ^ x ii x has rotation 
vector u). 

In the other direction, suppose the action is not free. Then there is a (k, I) ^ Jp^q and 
an a; £ Xp,q n X„ with Tk,ix = x. Clearly, such (fc, /) must satisfy {oj, k) + I = 0, that is 
{k,l) € lu,- Thus, Jp,g =^ lu,. □ 

The case that a; £ Q** is especially nice. We have the following: 

Proposition 3.11. cj G Q'' if and only if I^j has rank d. When (pi, gi), .... (pd, qd) £ I^j 
are linearly independent, thenXp^q is finite- dimensional andXp^q C X^^. In particular, when 
(pi, gi), . . . , (pd, qd) are principal periods, then Xp,g = X;^. 
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Proof. Let us suppose that uj G Q'', for instance uj = . . . , ^) for integers aj and bj. 
Then (a;, Pj) + Qj = for Pj := (0, . . . , 0, &j , 0, . . . , 0) and Qj = -aj. This shows that 
has rank d. 

On the other hand, when lui has rank d, then we can choose hnearly independent 
(pi, qi), . . . , {pd, Qd) & lui- If we now denote by p the d x d-matrix with integer coefficients 
(pi, . . . ,pd) and by g = (gi, . . . , qd) G Z'* the integer vector of length d, then we can write 
the equations {pj,uj) + qj — as the matrix equahty p'^ui + g = 0, where p^ denotes the 
transpose of the matrix p. This imphes that the final column in the rank-d matrix {p^ , q) is 
degenerate, i.e. that p^ is invertible. In particular, uj = —p^^q G Q'', with p~"^ the inverse 
transpose of the matrix p. 

Moreover, the fact that p is invertible implies that Xp,, is finite-dimensional. More pre- 
cisely, let us define 

Bp := p([0, 1)'') n Z"* . 

Then Bp is a fundamental domain for p, that is for every i G Z'* there is a unique k G Bp 
with k = i mod p(Z'*). It is not hard to show that this implies that the map s i-)- xja^ from 
Xp,q to is an isomorphism. 

Thus, we have that dimR'^'' = \Bp\ — vold(p[0, l)**) = | detp | and hence Xp,q is finite- 
dimensional. 

In turn this implies that Xp_q C X^^, because any x G Xp,, satisfies sup^^^d{xi — {uj, i)} = 
suPigBp{^i — ('^1*)} < oo and the configuration i i— >■ {ijJ,i} has rotation vector uj. 

If (pi, qi), . . . , (pd, qd) are principal periods, then the above implies that X^^ = Xt^nXp,, = 
Xp,^. □ 

After these general considerations, let us now return to the Birkhoff configurations defined in 
Section [3.1 1 Let us denote the set of maximally periodic Birkhoff configurations of rotation 
vector UJ by 

Bu := {x <E Buj \ Tk,ix = X when {k,uj) + I — Q } . 

The following theorem expresses that periodic Birkhoff configurations are automatically 
maximally periodic. 

Theorem 3.12. Let uj G R'*, denote c := lankz (Jt^) and let {pi,qi), . . . ,{pc,qc) G 1^ be 
linearly independent. Then 

Xp,, n Z3„ = B„ . 

Proof. Let x G Xp,^ DBui, that is x is Birkhoff, has rotation vector uj and Tp-,q-x = x for all 
J = 1, . . . , c. We need to show that whenever {uj, k) + 1 = 0, then Tk.ix — x. So let us assume 
that Tkjx 7^ X. Because x is Birkhoff, we may assume that Tk.ix > x: the case Tk^ix < x is 
similar. This assumption implies that T„k,nix — tJ^iX > x as well, for every n > 1. We claim 
that this is not possible. 

To prove this claim, we remark that there must be an n G N and mi, . . . , rric G Z so that 
n{k, I) = mj{pj,qj). This is because by assumption the {pj,qj) span a sublattice of lu, 
of maximal rank. We therefore have that Tnk,nix — (tj^^^^i o . . . o T^'^g^)x — x. This is a 
contradiction and hence, t^jx = x. □ 

In dimension d = 1, Theorem |3.12| simply says that a Birkhoff configuration of period 
{np,nq) automatically has period (p,q). That is, the period (p, (?) of a one-dimensional 
Birkhoff configuration can be chosen relatively prime. Theorem |3.12| is the d-dimensional 
variant of this statement. 



In spite of Theorem 3.12 it should be remarked that in general, B^i 7^ B,^, that is not all 



Birkhoff configurations of rotation vector uj are periodic. Counterexamples are easy to find. 



4 Classical Aubry-Mather theory 

We are now ready to discuss the most well-known results of classical Aubry-Mather theory in 
the context of lattice equations. The concepts and results of this section are widely known, 
but we chose to present them in a perhaps slightly unconventional manner. 



9 



4.1 Fully periodic minimizers 



Throughout Section 4.1 we will assume that u) G Q and (pi, q\), . . . , (pd, qd) £ Z x Z are 
linearly independent. 

We are interested in solutions to (2.8 1 that lie in Xp.g. We start by noting the presence 
of a variational structure. Recalling the definition Bp — U^r\ p([0, 1)''), one has 



Proposition 4.1. A configuration x £ Xp,q solves (2.8) if and only tf it is a stationary 
point of the periodic action function 

Wp.q : Xp,, ^ R defined by Wp,g{x) ■- Wb^{x) =Y1 ' (^^-^O) 

with respect to variations in Xp^g. 

Proof. We start by recalhng the shift-invariance of the local potentials, condition B, which 
says that Sj+k{x) = Sj{Tk,ix) for all k and / and all x € . Differentiation of this identity 
with respect to Xi then gives that diSj+k{x) = di-kSj{Tk^ix). These equalities respectively 
imply that for x £ Xp,, it holds that Sj+pk{x) = Sj{x) and diSj+pk{x) — di-pkSj{x) for all 
k G Z"^. 

Now let X G Xp^q, choose an i G Z'* and define G Xp,o by letting {ei)j = 1 if j = i 
mod p{Z'') and {ei)j — otherwise. Then x + ct € Xp,, and 



de 



Wp,q{x + eei) = J2 9^+pkSj{x) = X] 5Z 9iSj-pk{x) = J2 diSj{x). 

^0 kez''ieBp kez''jsBp j^zd 



Of course all these sums are finite. □ 

Note that Wp^q = is actually well defined for any a; G R^ , but in this section, we 

restrict it to a function on Xp,,. As such, it is shift-invariant: 

Lemma 4.2 (Shift-invariance). The function Wp^q is r-mvariant: for x G Xp,q and (k,l) G 
Z'* X Z arbitrary, it holds that Wp.q{Tk,ix) — Wp.q{x). 

Proof. In general, Sj+k{x) = Sj(Tk,ix), so that if x G Xp,,, then Sj+pk{x) — Sj{x) for 
all k G Z^. Thus, Wp,q{rk,ix) = E'.gB^ S,(rfc,a;) = E.gsl = E,sfe+SpS,(x) = 

EjgB ^j{x) ~ ^p,q{x)- The fourth equality follows as both Bp and k + Bp are fundamental 
domains of Z''/p(Z'*), so that for every j £ k + Bp there is a unique i £ Bp for which i = j 
modp(Z''). □ 



Theorem 4.3 (Existence). The action Wp^q attains its minimum on '. 



Proof. Since Wp^q{ro,ix) = Wp^q{x) for x G Xp,,, clearly Wp^q descends to a function on 
Xp_q/Z. Every element in this quotient space has a representative x with xg G [0, 1]. 

Choose a cube Cjv = {i G Z'' | < A'' for all fc = 1, . . . d} that contains Bp and choose 
a G Z'' and an n G N such that k + nBp in turn contains Cjv. Moreover, remember that 
Wk+nBp = n'^Wp^q on Xp,q. 

The coercivity of the Si, condition C, implies that for all j with ||j|| = 1, it holds that for 
X G Xp,q with xo G [0, 1], we have that limi^.^. Wk+nBp{x) = oo. And hence by induction 
that for all j G Cat and x G Xp^, with xq G [0, 1], it holds that limi^^^^ac Wk+nB^ix) = oo. 
Because Bp C Cn and Wp,q = n~'^Wk+nBp, this means in particular that for all j G Bp 
and X G Xp,q with xa G [0, 1] it holds that limi^,^ |_+cx) W^p,9(a;) = oo. Hence, Wp,q attains its 
minimum on Xp,,. □ 

The configurations that minimize Wp,q on Xp,, will be called p, g-minimizers. Note that 
other extremal points of Wp,q in Xp,,, such as saddle points, may also exist. Under certain 
mild conditions their existence will be proved later in this paper. 
The following lemma is well-known. We took the proof from |5]. 

Lemma 4.4 (Minimum - maximum property). Assume the periodic configurations x,y G 
Xp,g arep,q-minimizers. Then alsom := mm{x,y} andM := ma,x{x,y} are p, q-minimizers. 
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Proof. It is obvious that m,M G Xp,,. Write a :— M — x and P :— m — x and observe that 
a > 0, /3 < 0, while supp(Q) n supp(/3) = and y^M-^m — x = a + m = a + l3-[-x. The 
proof is done, if we show that 

Wp,,{x) + Wp,g{y) > Wp,g{M) + Wp,g{m). 

This is the same as showing 

Wp,g{x + a + l3)- Wp,,{x + a)- Wp,,{x + f3) + Wp,,{x) > 0. 

The left hand-side of this inequality can be put in integral form as 



^2 / /■! /■! 

di,kSj{x + ta + sl3)dsdt ) a-ifik 



Since supp(a) n supp(/3) = 0, we have that Qi/3i — for all i. Moreover, the twist condition 
di^kSj < for all i 7^ and the inequalities aiPk < 0, guarantee that the remaining terms 
in the sum are nonnegative. □ 

This is now used to prove the following famous lemma: 

Lemma 4.5 (Aubry's lemma). Assume the configurations x ^ y £ Xp,, are p, q-minimizers. 
Then either x <^y or y <^ x. 

Proof. We pursue a proof by contradiction. Denote again m := min{j:, y}. Suppose that for 
instance that m < x but that is not true that m <^ x. The case that m < y and not m <^ y 
is similar. The assumption implies that there are indices i, fc £ Z'* with ||i — fc|| = 1 such 
that mi = Xi and mt < Xk- Now we compute 



E f / di.iSj{tx + {l-t)m)dt^{xi-mi) . 



Recall that Xi = m^, while, by the twist condition, for every / 7^ i, it holds that di^iSj < 
and {xi — mi) > 0. Thus, every term in the above sum is nonpositive. But for the k chosen 
above, di^kSi < 0, while Xk — mk > 0. This proves that "^^j diSj{x) / diSj{m). This 
contradicts the fact that, by the lemma above, both m and x are p, g-minimizers and must 
therefore both be stationary. □ 

Corollary 4.6. Periodic minimizers are Birkhoff configurations. 

Proof. Let x £ Xp_q be a minimizer. Then for any fc G Z'' and / £ Z, we have that 
Wp^q{x) = Wp^q{Tk,ix) by the invariance property of Wp,q. This shows that also Tk,ix is a 
minimizer, whence, by the previous corollary, either Tk^ix <C x, Tk^ix — x or Tk^ix 3> x. In 
particular, a; is a Birkhoff configuration. □ 

Lemma 4.7. Let n £ N. Every p, q-mmimizer is an np, nq-minimizer and vice versa. 

Proof. Assume that x is an np, nq-minimizer, that is a minimizer of W„p,„q on X„p,„q. 
Then, by Aubry's lemma, x £ B„p^nq. Theorem |3.12| now implies that B„p,nq = Bp^q, so 
actually x £ Sp,, C Xp,,. Note now that on Xp,, it holds that Wnp,nq ~ n'^Wp^q and let 
y £ Xp,g C X„p,„,. Then Wp,q{x) = n~''Wnp,nq{x) < n~'^Wnp,uq{y) = Wp,q{y). Thus, a; is a 
p, g-minimizer. 

In the other direction, if a; is a p, g-minimizer and y is an np, ng-minimizer, then y £ Bp,q 
and W„p^nq{x) = n'''Wp^g{x) < n''Wp,g(j/) — W„p,„q{y), that is x is an np, ng-minimizer. □ 

The following result shows that p, g-minimizers are global minimizers. Recall that x is called 
a global minimizer if for every finite set B G Z'^ and every y with support in its r-interior 

one has that Wb{x + y) > Wb{x), with Wb{x) := Ejes ^ji^)- 
Theorem 4.8. Periodic minimizers are global minimizers. 
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Proof. Let x G Xp,g be a p, q-minimizer. If x is not a global minimizer, then there exists a 
finite set B C Z,'^ and a configuration y with supp(j/) C B^^\ such that Wb(2; + 3/) < Wb{x). 
Since B is finite, there exist & k £jf and an n £ N such that supp(j/) C -B C fc + -Bnp- Now 
define y G X„p,„g by setting t/i — yj when j is the unique point in + B„p for which j = i 
mod np(Z'^). In other words, y is the np-periodic extension of y\k+B„p- Then we conclude 
that 

Wnp,uq{x + y) - W„p,nq{x) = Wk+„Bp{x + y) - Wk+„Bpix) ^WB{x + y) - Wb{x) < 0, 

so x is not Tip, ng-minimizer. This contradicts Lemma [4.7[ □ 

Perhaps surprisingly, to prove the converse one needs to be slightly more ingenious. We have 
not found this statement anywhere in the literature: 

Theorem 4.9. If x £ Xp., is a global mimmizer, then it is a p, q-minimizer. 

Proof. Suppose that x £ Xp_g is not a p, g-minimizer. We will prove that this implies that 
X is not a global minimizer. Our assumption means that there is a i/ £ Xp,q for which 
< e := Wp,q{x) — Wp,q{y). This in turn implies that W„p,nq{x) — Wnp,nq{y) = n''e. 

By periodicity, we may assume that x <^y. Let us now define, for n £ N, the configura- 
tions X < y" < y hy 




yi if i G Bnp 
Xi otherwise 



Here, Bnp is the r-interior of Bnp. By definition, is a variation of x with support in this 
r-interor. It now holds that 

-ne+Y. {Sj{y)-SAyn) ■ 

Because the support oiy — y" is contained in Z''\BnJ and the range of interaction of the Sj is 
equal to r, the number of nonzero terms in the above sum is at most {2rf'^+'^\dBnp\ < En'^~^ , 
where _E is a constant depending only on r, d and p. 

Moreover, by compactness of [x,y\ ■= {z \ x < z <y}, there is a constant e > so that 
\Sj{y)\, \S]{y")\ < e. This then implies that 

WB„p(a;) - W^s„p(y") > en'' ~ 2Ee ■ t/-^ . 
Choosing n large enough, we see that x is not a global minimizer. □ 



4.2 Nonperiodic minimizers 

In this section, we show that global minimizers of all rotation vectors exist. They are 
constructed as limits of periodic minimizers. Moreover, we show that they satisfy a certain 
pairwise regularity. The results in this section are standard. 

Lemma 4.10. The set of global rmnimizers is closed in the topology of pomtwise convergence. 

Proof. Assume that x" is a sequence of global minimizers converging pointwise to 2:°°. Let 
i? C Z'* be a finite set and y a configuration with support in B^^K Then 

Wb^x" +y)-WB{x'')>0 . (4.11) 

But Wb is continuous with respect to pointwise convergence, so that taking the limit for 
n — ^ CO of equation (4.11 1, we find that Wb{x°^' + y) — Wfl(a;°°) > 0. So x°° is a global 
mirumizer. □ 

Theorem 4.11. For all rotation vectors cj G R'' and all local potentials Sj, there exists a 
global minimizer x G Bu) . 
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Proof. For any lj G R'*, we can take a sequence uj„ £ Q"^, such that hm„^oc i^n ~ oj, while 
(wn, k) +1 = for all the k and I for which {u),k) + 1 — 0. We take a corr esp ondi ng sequence 



(Pn, (Jn) of principal periods for which w„ := —p^^qn- By Theorems 4.3 and 4.8 there exists 
a global minimizer x" £ Bp^^q^ — Bui„ ■ In particular, x„ has rotation vector cun and satisfies 
Tk.ix" = x" for all fc and I for whic h (cj , k) + I = 0. Because the ljJ„ and uj lie in some 
compact subset K of R'', Proposition |3.6[ guarantees that there is a subsequence of the x" 
that converges pointwise to a Birkhoff configuration x°° C Bk- By continuity of the rotation 
vector X i— >■ uj{x), see Proposition |33j 2;°° actually has rotation vector w. Moreover, the limit 
x°° will have the same periodicities: denoting the converging subsequence also by x", the 
continuity of rk,i implies that Tk,ix°° = rfc,i(lim„^oo a;") = hm„^oo Tfc,;a;" = lim„-^oo a;" = 



x°° for all k and I with (a;, k) + 1 = Q. Finally, is a global minimizer by Theorem 4.8 and 



Lemma l4T0l □ 

The following result expresses the regularity of pairwise comparable stationary solutions. It 
is the analogue of a Harnack inequality for elliptic PDEs. 

Theorem 4.12 (Elliptic Harnack inequality). Let x < y be two Birkhoff configurations with 
rotation vector in the compact set K cM.''. Suppose that x and y are stationary for the local 
potentials Sj. Then there is a constant S, depending only on K and ||j — fc||, such that for 
all i and k, 

iUk - Xk) < S{y^ - Xi) . 
In particular, if x < y, then x <^y. 

Proof. By interpolation; let x and y be stationary and Birkhoff and let i,k and assume 
first that ||j — = 1. Choose a B with i £ B''^-' and recall the definition Wb{x) = 
"^jeB ^ii^)- Then, by stationarity. 



=d,WB{y) - ^^WB{x) = f ^\ 9,Sj{Ty + (1 - r)x) dr = 

E (fd,,iSAry + {l-r)x)dT\ {yi^xi) . 

Since, by the twist condition C, the only possibly positive terms on the right hand side are 
the {di^iSj{Ty + (1 — t)x)) (yi — Xi), the right hand side is less than or equal to 



di^tSj{Ty+(l-T)x)dTUyi~Xi)+i ^ di,kSi{Ty + {l-T)x)dT\ (yk-Xk) 



Now, because x and y are Birkhoff, so is every ry + {1 ~ t)x and hence by Corollary 3.7 
there are constants A, C > 0, depending only on the compact set K, such that for all j and 
all ||i — fc|| = 1, it holds that di^kSi < —A, while di^iSj < C for all i and j. Thus, 

< {2rfC{y^ - X,) - 2d\{yk - Xk) . 

This proves the theorem for ||i - fc|| = 1 with 5 ^ 5i := {2rYC/2dX. For ||z - k\\ > 1, the 
result then follows by induction and it holds for 5 = □ 

4.3 Aubry-Mather sets 

We make the following definition: 

Definition 4.13. An Aubry-Mather set Ai C R^ is a collection of configurations with the 
following properties 

• is nonempty and closed under pointwise convergence 

• AI is strictly ordered, i.e. for every x,y(^M,x<^y,x — yorx^y 

• AI is shift-invariant: if 2: £ AI, then for every (fc, £ Z'* x Z, also rk,ix £ AI 

• Every t £ AI is a global minimizer of the variational recurrence relation ( |2.8[ ) 

• A4 does not contain any strictly smaller set with the properties listed above 
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The strict ordering and the shift-invariance of an Aubry-Mather set A4 imply that any 
configuration x £ A4 is Birklioff and hence has a rotation vector cj — ijj{x). The ordering of 
A4 moreover implies that this rotation vector is independent of the choice of x £ M, that is 
oj — uj{M) and thus, M C Bi^. 

Recall that Theorem |4. 11| states that for every rotation vector uj there exists a minimizer 
X £ Bui for which Tk,ix = a; as soon as {Ld,k) +1 — 0. This in fact implies that a certain 
Aubry-Mather set A4(x) C Bu, exists. This A4{x) is constructed as follows. One starts by 
defining the collection M{x) C Bui as the closure with respect to pointwise convergence of 
the set of translates of x: 



M{x) ■- {Tk,ix j {k,l) G Z'i X Z }. 
This is almost an Aubry-Mather set: 

Lemma 4.14. Let x G Buj be an action-minimizer with the property that Tk^ix — x when 
{ui, k) + 1 = 0. Then M{x) ts nonempty, closed, strictly ordered, shift-invariant and consists 
of minimizers. Moreover, for every y £ M (x) it holds that r^^iy — y as soon as {uj, k}+l = 0. 
When Lu £ Q'*, then M{x) is an Aubry-Mather set. 

Proof. By definition, Mix) is nonempty and closed. 

We note that when x is a minimizer, then so is Tk,ix and because any pointwise limit of 



minimizers is a minimizer itself, by Lemma 4.10 we see that M(x) consists of minimizers 
only. _ 

When y £ JVl[x), say y = Vmin^ooTk^^^i^x, then the continuity of t^.i implies that 
Tk,iy = Tfe,;(Um„^oo Tfc„,;„a;) = lim„^oo rk,i{Tk„,i„^x) = lim„^oo Tk+k„^,i+i^x and thus, M{x) 
is shift-invariant. 

The fact that a; is a Birkhoff configuration means that the collection {rk^ix \ {k,l) £ 
Z'* X Z } is ordered. Now let y and z be elements of M{x), say y = Ymin^aa '''k„,i„x and 
Z = limn^oo TK,^,L„x. We claim that y < z or z < y. If not, then there are i,k with 
j/i < Zi and yk > Zk. The pointwise convergence then implies that there are n and A'^ so that 
< iTK^.L,^x )t and {Tk„,i„x)k > iTK„,LMx)k. This is a contradiction. The second 
conclusion of Theorem 4.12 now implies that y <^ z, y = z or y ^ z, that is A4{x) is strictly 
ordered. 

The penultimate conclusion of the lemma follows from the continuity of Tk^i and the 
fact that Tk,ix = x when {uj,k) -\-l = Q. Namely, for such k and I and for y £ M(x), say 
y = lim„^oo T-fc^.i^i, we have that Tk,iy = rfc,;(lim„^oo Tfc„,i„a;) = lim„^oo 'rfc,i(Tfe„,;„a;) = 
lim„^tx, rk,^,i,^{Tk,ix) = hm„.^oo Tk„,i„x = y. 

Finally, when cj £ Q'', then our assumptions imply that x is periodic, say x £ Xp,q, with 
(p, q) a collection of principal periods for uj. This implies that the r-orbit of x is finite. Thus, 
A4{x) is equal to this single r-orbit and cannot contain any proper nonempty r-invariant 
subset. □ 



It is clear from the proof of Lemma 4.14 that when lo £ Q'', then every Aubry-Mather set 
is finite and consists of the translates of one periodic minimizer. Thus, the Aubry-Mather 
sets of rational rotation vector do not need to be unique. 

On the other hand, when UJ £ R'^\Q'^ is irrational, then Ml (x) may fail to be an Aubry- 
Mather set. Then one replaces M{x) by its recurrent subset 

M{x) := {y £ M{x) \y = lim Tk^^i^^y for a sequence (fc„,/„) with {uj,kn) + In 7^ } . 

Before proving that this M{x) is indeed an Aubry-Mather set, let us define for a configuration 
y £ M{x), the configurations 

y~ ;= snp{Tk.iy < y} and j/+ ~ inf {Tk,iy > y} . 

We remark that, by definition, y £ M(x) if and only if y = j/~ or y = i/^, or both. We now 
have the following technical result: 

Proposition 4.15. For y,z £ M{x) it holds that y~ = sup{Tk,iz ^ y~} and y"*" = 
mi{Tk,iz > y+}. 
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Proof. Let us prove the first equality: tiie proof of tlie second one is similar. We denote 
z~(y^) := sup{Tfc,;z ^ y~} and we argue by contradiction. That is, we suppose that 
z~(y~) 7^ , and hence, that z~{y) <^ y~ . Then, because y~ can be approximated from 
below by translates of y by definition, there are k and I so that z~ (y~) <^ Tk,iy <C y - This 



implies that r^^iy <^ y and in view of Proposition 3.8 we must therefore have that (co, k)+l < 
0. Applying T-k.-i to the inequality z~ {y~) -ti Tk,iy, we obtain that T-k,-iz~ {y~) <C y. But 
because (cj, — fc) — Z > 0, we must also have that z~{y~) <C T-k,-iz~ (y~). Hence, z~{y~) <C 
T-k,-iz~ (y^) <^ y. But this contradicts the definition of z~{y~), because by continuity of 
T-k,-i, if z~{y~) = hm„-+cx) Tfc^j^z, then also T-k,-iz^ (y^) = hm„^oo T_fe+fc„,_i+i„z is a 
limit of translates of z that lie below y~ . □ 

We are now ready to prove: 

Theorem 4.16. When cj G R'^/Q'', then M{x) is the unique Aubry-Mather set contained 
in M(x). 



Proof. Proposition 4.15 says that any y G A4{x) is a limit point of the r-orbit of any 
z £ A4{x). Thus, any nonempty, shift-invariant closed subset of A4{x) should contain 
A4{x). It remains to show t hat A4 (x) is nonempty, shift-invariant and closed. 



First of all. Proposition 4.15 applied to z = y~ and z = y respectively, says that 
{y~)~ = y~ and (y"*")^ = y'^ , i.e. that y~ and y"*" are recurrent. This shows that A4{x) is 
nonempty. 

Shift -invariance of Adi^x^ follows from the continuity of t^.i'. when y — limj^— >oo '^kn.iny-; 
then Tk,iy = lim„_j.oo 'rk„,i„{'rk,iy)- 

To prove that M{x) is closed, assume that lim„^oo J/n = y pointwise for a sequence 
j/n of recurrent configurations. When the limit y is not recurrent, then y~ ^ y ^ y"*", so 
that there is an n for which y~ <gi y„ <Si y^ . But y„ is recurrent, hence y„ 7^ y, while by 
Proposition |4. 15] y„ can be approximated by translates of y. Hence, there are k and / such 
that y~ <^ Tk,iy ^ y'^ and Tk.iy 7^ y. This contradicts the definition of y~ or y"*". □ 

Remark 4.17. A theorem of Bangert 2, in the case of elliptic PDEs, states that when 
uj £ R'*\Q'', then the recurrent subset actually does not depend on the choice of the Birkhoff 
minimizer x £ Bui. In other words, that when x,y £ Bu are such that Tk.ix = x and Tk,iy = y 
whenever (ui, k) + I = 0, then M{x) = M{y). 

The proof of this theorem is nontrivial. The essence of it lies in proving an Aubry lemma 
for recurrent minimizers, that is to show that if i £ A4{x) and y £ A4{y) are recurrent, then 
X <^ y, x = y or x ^ y. 

We claim that a similar theorem holds for lattices instead of PDEs, but we will not prove 
this, as it is not essential for the remainder of this paper. As a result, the Aubry-Mather set 
of an irrational rotation vector is unique. 

The following well-known result shows that the set of recurrent minimizers can have a 
complicated topology. We recall that a topological space C is called a Cantor set if it is 
closed, perfect and totally disconnected. "Perfect" means that every element c £ C is a limit 
of points in C\{c}. "Totally disconnected" means that for any two elements ci,C2 £ C one 
can decompose C as the disjoint union of closed sets Ci and C2 with ci £ Ci and C2 £ C2. 

Theorem 4.18. If u) & R''\Q'*, then M(x) is either connected or a Cantor set. 

Proof. The recurrent subset is perfect by definition: for every y £ A4{x), it holds that 
y = limn-^oo '^k^,i^y, where by Proposition |3.8| the condition that (cj, kn) + ln / guarantees 
that Tk^,i^y 7^ y for all n. 

We will now show that when Mix) is not connected, then there is a j/ £ Mix) so that 
y~ 7^ 1/^. So let's assume that Mix) is not connected and write Mix) = U VJV for two 
nonempty closed subsets U and V with U C[V — ^. We may assume that there exist u £U 
and V £ V so that u <gi v, whence we can define y :— sup{M £ 17 jit <C v}. Clearly, y £ U, 
because U is closed. Hence, y <^v. We claim that y^ 7^ y. This is easily proved: if y"*" = y, 
then y — lim„_>oo n^^i^y for a sequence with y ^ '''k„,i„y «C v- By definition of y, it must 
hold that Tk„,i„y £ V. Hence, because V is closed, also y £ V, which is a contradiction. 

The next step is to observe that an order interval [y~ , y'^] :— {z £ M.^ I J/~ < < y^} 
can never contain any recurrent elements other than y~ and 3/+. Namely, if y~ <^ v <^ y^ 
were such a recurrent element, then by Proposition |4. 15[ it can be approximated by translates 
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of y, so that there are k and I with y~ <^ Tk,iy ^ y'^ ■ This contradicts the definition of y~ 
or y^ . This is why we call the order interval [y~ ,y^] a gap in the Aubry-Mather set. 

Now we show that when A4{x) is not connected, and hence contains at least one gap 
[y~ ■,y^]j then between any two elements w,z £ A4{x) there exists a gap. Namely, for any 
given pair w <^ z, either [w,z] is a gap, or there is a recurrent element w <^ u <^ z. 



By Proposition 4.15 this u can then be approximated by the r-orbit of y~ , which implies 
that there are k and I so that w <C Tk,iy~ <C z. But when [y~,y^] is a gap, then so is 
[Tk,iy~ ,Tk,iy^]; since Tkj is order-preserving. We must therefore have that w ^ T^^y" ^ 
n.iy^ < z, i.e. that there is a gap between w and z. 

This implies that M{x) is totally disconnected: if w,z £ M{x) with w <^ z, then there 
is a gap [y~,y^] with w < y~ <^ y'^ < z and hence M{x) splits as the disjoint union of 
the closed sets {u G A4{x) \ u < y~} and {v £ M{x) \ v > y^} that contain w and z 
respectively. □ 



The proof of Theorem 4.18 shows that for any y G A4{x), in the order interval 

~{zeR^' \y' <z<y+} 

only the elements y~ and y^ are recurrent. Hence, when y~ 7^ y"*", then [j/",?/"''] is called 
a gap in the Aubry-Mather set. Moreover, in the case that M{x) is not connected, then 
between any two recurrent configurations there exists such a gap. 

When M{x) is connected, then we say that it forms a foliation: for every i £ 11'' and 
every ^ € K there is a unique y £ Mix) so that yi = ^. In the case that Mix) is a Cantor 
set, one says that it forms a lamination: for every i and every ^ there is at most one y so 
that yi = 

Both foliations and laminations by minimizers occur in examples, for instance that of the 



Frenkel-Kontorova lattice ( 1.1 1. In fact, when V{£_) = 0, then the Aubry-Mather sets are all 



of the form Ai{x'^' ) := {a;"'^ | ^ G R}, where we recall that the linear configuration 2;"'^ is 
defined by x"'* = {uj, i) + These Aubry-Mather sets are clearly connected. 

On the other hand, the following theorem says that when the onsite potentials V^(^) are 
sufficiently oscillatory, then the Aubry-Mather sets must be Cantor sets: 

Theorem 4.19. Let Sj be local potentials satisfying conditions A-E and let K C M.'' be 
a compact set. Then there exists a number M > 0, depending on the Sj and on K, 
such that for every 1-penodic twice continuously differentiable function V — V{^) with 
osc V := max5_^gii(V(^) — Vii^)) > M, the collection of local potentials Sj defined by 
Sj{x) = S j(x) -\-V (x j) does not possess any connected, strictly ordered shift-invariant family 
of global minimizers of rotation vector uj £ K . 

Proof. Because Bk/1 is compact and the functions Sj are r-invariant and continuous, their 
oscillation over Bk is bounded and, say, equal to := osc Bk^j ~ '^^'^x.yeBKi^ii^) ~ 
Sj{y)). Let M > (2r -I- I)'* A'', where r > 1 is the finite interaction range of the local 
potentials Sj , and choose a smooth 1-periodic onsite potential V with oscillation larger than 
M. Assume for instance that V{^) — V{v) > M for certain ^, £ R. 

We will now prove that if a configuration x £ Bk has xo = then it can not be a global 
minimizer. In other words, that a; is a "gap configuration". This is easily shown by defining 
J/ : Z'* — >■ R by setting = for i 7^ and 2/0 = 1^ — ^. Now choose a finite subset i? C Z'' 
such that Bq C B. Then supp y — {0} C -B^''' and we compute that 

wb{x) ~ wb{x + y)^Y. ^^■(^) - ^^(^ + 2/) = E •^^■(^) - *^^(^ + y) = 

jeB j&Bi 
- V{v) + Sj{x + y) - Sj{x) > M - (2r + 1)''N > 0. 

This shows that x is not a global minimizer. □ 
Example 4.20. For the Frenkel-Kontorova lattice. Theorem |4.19| can be improved upon 



considerably. In fact, by Lemma 3.5 osc Bu, i^j ^ Xk) < 2, which is independent of uj. There- 



fore, the oscillation over B of the interaction potential ^ X]||j_fcji=i(^j ~ a^fe)^ is bounded 
above by 1. Thus, for any onsite potential V{S,) with oscillation larger than 2d, the Frenkel- 
Kontorova lattice with local potentials Sj(x) = ^^^\k-j\\=ii^3 ~ a^fc)^ + Vi^j) does not 
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have a connected family of global minimizers of any rotation vector at all. 

The latter result for the Frenkel-Kontorova lattice is well-known in dimension d — 1. It 
turns out that the one-dimensional Frenkel-Kontorova lattice is equivalent to the Chirikov 
standard map Tv, see Appendix A. As such, Theorem |4.19| and the discussion above say 
that for any onsite potential V with oscillation larger than 2, the standard map Tv has no 
rotational invariant curves. In the case that V has the "standard" form V{£,) = cos(27r5), 
so that osc V = we obtain that there are no rotational invariant curves for k > Svr^. In 
fact, in this case the much stronger computer-proved bound fc > || is actually known, see 



5 A formal gradient flow 

The idea of studying globally stationary solutions by means of a formal gradient flow goes 
back to Gole, see [7] . We will review his ideas in this section. The new result is a parabolic 
Harnack inequality, see Theorem |6.4| 

The study of the formal gradient flow starts with the observation that one can assign a 
meaning to the partial derivatives of the formal, and generally nonconvergent sum W{x) = 
X^jGZ'i '^i(^)' namely as follows. Since the potentials Sj are of finite range, every variable Xi 
appears only in finitely many terms of the formal series. Hence, we may write, with a slight 
abuse of notation, 

(VVK(a;)), := d^W(x) = ^ ^^S,{x) . 

I I <»" 

Note that VW : -> K' is well-defined as soon as the Sj are continuously differentiable 
and that VW is the formal gradient of W with respect to the Z2-inner product {x,y)2 = 

We remark that x is globally stationary if and only if \/W{x) = 0. In this section, we 
shall nevertheless view such x as stationary points of the auxiliary differential equation 

f - • 

This differential equation shall be defined on an appropriate Banach subspace X C of 
configurations, for which its initial value problem has existence and uniqueness of solutions. 
The corresponding flow is called the negative gradient flow of W. The motivation to study 
the negative gradient flow is simply that it will help us flnd globally stationary solutions. 

The Banach subspace we choose to work with is the exponentially weighted conflguration 
space 

X := e R^'l 1|.t11x E < ~> ^ ' 

where we recall that := Ylt^i First of all, the space of Birkhoff conflgurations is 
contained in X; 

Lemma 5.1. S C X. 

Proof. This follows because every x £ B has a rotation vector, say lu, and \xi—xo — {oj, i)| < 1. 
This implies that \xi\ < ■ + \xi,)\ + 1 and hence 



□ 



We moreover note that the topology B inherits from X is exactly that of pointwise conver- 
gence: 

Proposition 5.2. Let a; G X and for all n G N, let x„ G X. Then lim„_j.oo \\xn — a;j|x = if 
and only if lim„_K3o Xn = x pointwise. In particular, a sequence in B converges in X if and 
only if it converges pointwise. 

Proof. The flrst claim is obvious. The second claim follows because S is a closed subset of 
X. □ 
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Before showing the existence of the negative gradient flow on X, we need the following simple 
lemma, which shows that the shift maps r^,; are Lipschitz on X: 

Lemma 5.3. Let x,y € "K and (k,l) £ Z'' x Z. Then Tk.ix £ X and rk,iy € X, while 
||Tfc,,a:-rfe,,y||x < 2ll*"||x-2/||x. 

Proof. First of all, 

II i,u iiA / y 2ll»ll / ^ 2ll*ll + ll''ll ~ ' 2ll'+*ll 

iijA iez'i iez.'' 

Therefore, ||rfc,ia;||x = | |rfc,o2: + ^| |x < ||Tfe,ox||x + |K||x < 2"'°" ||j:||x + ||i||x < oo and similarly 
ior Tk,iy. In particular, ||rfe,;a:: - rfcjy||x = ||rfc,()(a:: - y)||x < 2ll''ll Hx - y||x. □ 

In particular, this means that Tk,t : X — > X is continuous in the topology of pointwise 
convergence: if x" — >■ x°° pointwise, then Tk,ix" — > Tk,ix°° pointwise. Of course, this is also 
clear without Lemma 15.31 

The main result of this section is the following theorem, which says that under the 
condition that the local potentials Sj are twice continuously differentiable with uniformly 
bounded second derivatives, then —VW indeed defines a fiow on X. Moreover, this flow has 
the regularity properties one expects it to have. 

Theorem 5.4. Assume the local potentials Sj satisfy conditions A, B and E, that is they 
are twice continuously differentiable with uniformly bounded second derivatives, they depend 
on finitely many variables and are shift-invariant. Then the vector field —VW : X — >■ X is 
globally Lipschitz continuous, i.e. there is a constant L > 0, depending only on the constant 
C of condition E and the interaction range r of condition A, such that for all x,y £ H, 

\\VW{x)~VW{y)\W<L\\x~y\\^. 

Hence, the initial value problem ^ = —'VW{x), x{0) = xo on X has global-in-time existence 
and uniqueness of solutions and defines a complete flow 1 1— >■ on X. This flow is Lipschitz 
continuous, i.e. there are constants Lt > 0, depending only on L, such that for all t (I'M. and 
x,y €X, 

11**2;- *tyjix < Lt\\x-y\\x. 
Moreover, this flow depends Lipschitz continuously on —VW. This means that there are 
constants Lt > 0, depending only on L, such that for all t £ M. and a;, y G X and for all 
—VW and — VW with Lipschitz constants < L and respective complete flows 'i't and ^t, 

sup|!*ta;- *t2;|jx < sup ||Vl¥(a;) - VW'(a;)||x. 

xex lex 

Proof. Using the uniform bound that \di,kSj \ < C, see condition E, we will prove that —VW 
maps X to X and is globally Lipschitz continuous. The usual ODE theory then provides the 
existence of a complete flow 1 1— >■ ^tj on X. 
Thus, let a;, y e X. Then first of all 

\-VW(x),+VWiyU< J2 \^^S,{y)-^,SJix)\= ^ J' (^Q^s.iry + {1 - r)x)) dr 



\\j~i\\<r 



y di^kSjiry -\- {1 - T)x)dT 



\yk-Xk\< C" y y \^k-yk\ 

ll*— ill<'' \\j — i\\<r 



But this implies that 



||-VM/(:c) + VM/(y)||x<Cy 2-ll"ll y y \xk - yk\ ^ 

igZ<* I |fe-j| l<r| |j — 1| |<r 

C y y y 2~"'"|a;i+m+„ — = C y y ||rm+„,oa; — Tm+n,oJ 

||m||<r||n||<rigZ^ ||m||<r||n||<r 



By Lemma 5.3 and the fact that in the sum above | |m + n| | < 2r, we know that | |rm+,i.oa; — 
Tm+n, oyj|x < 2^'Hx — y||x. Heucc, noting that |{i G Z"* | ||i|| < r}| < (2r + 1)'*, we obtain 
that 

\\-\/W{x) + \7Wiy)\\x<L\\x-y\\x , 
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where L ■- 2'^''C{2r + 1)^'*. On the one hand, this shows that —VW is giobaUy Lipschitz 
continuous. On the other hand, choosing y = 0, we see that | \—\/W{x)+\/W{0) ||x < -LUsjjx, 
or II - VW(x)\\x < L\\x\\x + ||VW(0)||x, that is -VW maps X into X. 



This implies the existence and uniqueness of solutions of the initial value problem 



dx 



-\7W{x), x(0) = xo in X, that is the existence of flow maps *t : X X for all t & R. The 
Lipschitz continuity of ^'t follows, as usual, from an application of Gronwall's inequality: 
first one shows that \\'^tx — ^tj/Hx < — y\\x + L jj'' ||^'t2: — '^Ty\\xdr. This then implies 
that ||*ta; - *ty||x < Lt\\x - y\\x, with Lt = e-^l*L 

For the last part of the theorem, let VW and VW be two vector fields with Lipschitz 
constants < L and complete flows and respectively. Call x{t) = ^tX and x{t) = 'i'tx. 
We then have 

\\x{t) - i(t)||x < ||VW^(i(r)) - VW{x{T))\\^dr < \\\7W{x{T))-\7W{£{r))\\^dr+ 
Jo Jo 

\\\/W{x{t)) - \/W{x{T))\\xdT<\t\ sup ||VM/(a;)-VW^(x)||x + L /' ' ||i(r) - x(r)\\xdT. 

x£X Jo 

Thus, by Gronwall's inequality, ||*ta:: - *tx||x < Lt sup^g^ ||VW(a::) - VW'(a;)||x with Lt = 
|t|e^l*l. □ 

Remark 5.5. It is not true in general that —VW : X — >■ X is a map. Hence, contrary 
to a claim made in the ^'j in general can not be assumed either. 

By Proposition |5.2[ the first part of Theorem |5.4| implies that VW : X — >■ X is continuous with 
respect to pointwise convergence: if lim„_s.oo 2;" — x°° pointwise, then lim„_>oo \7W{x") = 
VW(a;°°) pointwise. 

Similarly, the second part of Theorem |5.4| implies that for every t £ R the flow map 
^t : X — >■ X is continuous with respect to pointwise convergence. 

Part three of Theorem |5.4| implies that if VW",VW°° : X — >• X is a sequence of for- 
mal gradient vector fields with a uniform Lipschitz constant and corresponding flow maps 
^T, : X X and such that \/W" VW°° uniformly on X, then for aU t also -f]^ ^ 
uniformly on X. 

In the remainder of this section, we will formulate a concept of convergence for a sequence 
of flnite range potentials S" that guarantees that their corresponding gradient vector flelds 
and flow maps converge uniformly. It turns out that it is enough to require the convergence 
of the gradients of the S". We will flrst of all need to define what it means for collections 
of gradients of finite range potentials to be "close" . Remembering the definition in Sec- 
tion |2] of the partial derivatives dji,...j^S of a fc times continuously differentiable function 
S : — >■ E of finitely many variables, we now define: 

Definition 5.6. Let S : — > R be an m+1 > times continuously differentiable function 
of finitely many variables, that is S(x) = s{x\b) for a certain finite subset B G Z'^ and an 
m + 1 times continuously differentiable function s : R^ R. Then we define the uniform 
C""(R^ ) norm || VS||^,„jjjz£ij G [0, 00) of the gradient of 5* as the finite sum of suprema 

l<fc<m-|-l jj,...,j^gZ'' xeR'^'^ 

We note that if Sj : R'^ — > R is a collection of m -I- 1 times continuously differentiable, 
shift invariant finite range potentials, that is if Sj{x) — Sj{x\B':) for some m + 1 times 

continuously differentiable function sj : R^j — >■ R and Sj{Tk,ix) = Sj+k{x) for all j,k and 
I, then l|VSi||p„^jjj;dj = llVSjU^^^jj^dj for aU i,j G Z"*. With this in mind, we first of afl 
prove: 

Proposition 5.T. Let Sj, 5j : R^ — )■ R be two m -f 1 > 1 times continuously differentiable 
and shift-invariant collections of finite range local potentials, say Sj{x) — Sj{x\B^) and 

Sj{x) = Sj{x\B':) o,nd denote their corresponding gradient vector fields by VW and VW . 
Then, there is a constant L > 0, depending only on the dimension d, such that 

sup||VVy(:r) - VW^(a;)||x < I||VSo - VSoH^o^z^,- 
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Proof. We have that 

|iVW^(x)-VH/(x)||x< ^2-ll'll Yl \^^S,{x) - d.S,{x)\. 
iez<^ l|j-i||<'- 

By shift invariance, diSj{x) — di-jSo{Tjflx) and similarly for Sj, so that sup^ \diSj{x) — 
diSj{x)\ — sup^ \di-jSo{x) — ^i-jS^^(x)\, and consequently 

supi|VW^(x)-VVK(x)|ix< ^2~ll'll ^ sup|9,_,So(x)-9._,So(2:)l = 
^2-ll'll ^ supi9,SoW-a,So(x)l <IllVS'o-V5o||po(jj,.), 

withi; = E>ez<*2'"'". □ 
We are now ready to define what it means for a sequence of local potentials to converge: 

Definition 5.8. Let S^,S^ : — >■ R be a sequence of collections of m + 1 > 1 times 
continuously differentiable, shift-invariant functions of finite range r. Then we say that the 
VS" converge to the VSf uniformly in C""(M^'') as n ^ oo if 

With this definition, we can then prove the following corollary of Theorem |5.4| It trivially 
follows from our definitions. Theorem |5 . 4| and Proposition |5.7[ 



Corollary 5.9. Let : — > R 6e a sequence of continuously differentiable local 

potentials of finite range r, with corresponding gradient vector fields VW" and VW°° , and 
assume that VS"" VS°° uniformly m C°(R^''). Then VW" ^ VW°° uniformly, i.e. 

lim sup||VW"(x)-VW°°(2:)|iy = 0. 

Moreover, m the case that the S" and are twice continuously differentiable with uniformly 
bounded second derivatives, so that — VW" and — VM^°° have well defined flow maps '5" and 
"^f, then it also holds for every t G R that — 5> '^f uniformly, i.e. 

lim sup I l^/^a; — ^tj^a;! |x = 0. 

n^co j.gx 

6 Properties of the gradient flow 

In this section, we collect some qualitative properties of the formal negative gradient flow 
that was introduced in the previous section. 

First of all, not surprisingly, it is equivariant with respect to shifts: 

Proposition 6.1. Let {k, I) e Z'' x Z and t £ R. Then o r^,; = r^j o 

Proof. By the shift-invariance of the local potentials Sj of condition B above, we have that 
Sj-k{Tk,ix) = Sj{x) for all k,l and j. Differentiating this identity with respect to Xi+t, we 
find that diSj-k{Tk,ix) — di+kSj{x). Assume now that ^ — — {\/W{x))^ for all i. Then, 

\\j-(i+k)\\<,- 
- ^^SJ-k{rk,lx)=-{VW{rk,lx))^ . 

\\(j-k) — i\\<r 

In other words, when t x(t) is a solution of the negative gradient fiow, then so is f 

Tk.ix{t). □ 



Proposition |6.1| implies in particular that the spaces Hp.q of periodic configurations are 
invariant under the gradient flow. 

The following well-known property of the negative gradient flow is the analogue of the 
comparison principle for parabolic PDEs, cf. [7] or [TT]. It is a direct consequence of the 
monotonicity condition D. 
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Theorem 6.2 (Strict monotonicity of the parabolic flow). Let a;, y G X such that x < y. 
Denote by the time-t flow of x — —\/W(x). Then for every t > 0, "i/tx <C ^'tj/- 

Proof Denote x{t) = "i/tx and y{t) = and define u(t) := y{t)-x{t). Note that u(0) > 
and that u satisfies the following linear ODE: 



u^{t) = -^M{y{t)) + ^M{x{t))^ — Y. d,SATx{t) + {l^T)y{t))\dr^ 
[f -^^,kSJ{rx(t) + (1 - T)y{t))dT] Uk{t) - {H{t)u{t)), . 

Here, for every t, the operator H{t) is Lipschitz from X to X, by a proof similar to that of 
Theorem |5.4[ Recall that di^kSj < when i ^ k, whereas di^iSj < C. This implies that 
there is a constant M > such that the operators H{t) := H{t) + Mid : R^'' R^'' are 
positive: u > implies H{t)u > 0. 

Note moreover that both the H{t) and the H{t) are uniformly bounded operators, whence 
the ODEs ii = H{t)u and v = H{t)v define well-posed initial value problems. More impor- 
tantly, u{t) solves u = H{t)u if and only if v{t) :— e^''^u{t) solves ii = H{t)v. We will now 
prove that for every t > and every i, Vi{t) > 0. Then, obviously, Ui{t) > as well, which 
then proves the theorem. 

To prove the claim on v{t), we solve the initial value problem for v — H{t)v by Picard 
iteration, that is we write 

v{t)=(Y,H^-Ht)\v{0), (6.12) 

\n=0 / 

where the H^"^{t) are defined inductively by 

H^-^Ht) = Id and H^"\t) := /' H{i) o //(""^'(t) di forn > 1 . 



Observe that the positivity of H{t) implies that the H^"\t) are positive as well. Because 
11(0) = m(0) > 0, we can therefore estimate, for any i, fc £ Z'' with ||i — — 1, 

v,{t) = (t)^;(0)) > {H{t)vm^ > 

\n=0 / i 

{Jo Jo + " T)y{i))dTddj vu{Q) . (6.13) 

Now choose a. k £ such that Ufc(O) > and recall that di^kSi < 0. Then from (6.131 it 
follows that if I |i — fell = 1, then for aU t > 0, Vi{t) > 0. 

To generalize to the case that ||i — fe|| 7^ 1, let us choose a sequence of lattice points 
fe = io,...,ijv = i such that — in-i\\ ~ 1 and = II* ^ ^11- Then, by induction, 
V (^)^ > for aU n. Thus, if Vk{0) > and t > 0, then Vi{t) > 0. □ 



Theorem |6.2| immediately gives us the following important corollary. 

Corollary 6.3. Let to £ W^. Then Bu, is positively invariant under the negative gradient 
flow: "iltiBuj) C Bu,, for every t > 0. 

This just follows because the strict monotonicity of the parabolic fiow implies that 'i't pre- 
serves the inequalities that define B^. 

The following is a quantitative version of Theorem |6.2| It will be crucial in the remainder 
of this paper and we have not found it elsewhere in the literature. 

Theorem 6.4 (Parabolic Harnack inequality). Let t > Q, K <Z W'' a compact set and 
x,y £ Bk ■~ y-iuieKBuj such that x < y. Then there exists a constant L > 0, depending 
on K, \\i — k\\ and t, such that for all i,k € Z'*, 

(*ty)i - (*ta;)i > L{yk - Xk) ■ 
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Proof. The proof is a quantitative variant of the proof of Theorem |6.2| We start by recalhng 
that by Corollary |3.7[ there is a constant A > 0, depending only on K, such that di^kSi{z) < 



— X < for all ||i — fc|| = 1 and z G UuiskBi^j. Then (6.131 shows that if ||i — = 1, then 
Vi{t) > LiVk{0), with Li — tX. 

To generalize to the case that ||A: — i|| 7^ 1, we again choose a sequence of lattice points 
k — io, . . . , ijv ~ i such that \\i„ — in-i|| = 1 and A*' = ||i — k\\. Then there is a constant 

= ^ depending only on Jf, i and I |i — fc| I such that ?; (^). > L'f^^v ^ ^"^^''^ ^ for all 

n. Thus, v^{t) > LNVkiO) with Ln = (L'at)^- 

This proves that Ui{t) > Luk{0) with L = e-^^^Ln = e-'^''{Xt/\\i - /c||)ll'-'=ll. □ 



Note that for Birkhoff configurations, both the strict monotonicity, Theorem |6.2| and the 
elliptic Harnack inequality. Theorem |4.12| follow directly from this parabolic Harnack in- 
equality. 

We moreover remark that under the uniform twist condition that di^kSi{z) < —A < 
for all z G X and ||i — fc|| — 1, the above parabolic Harnack inequality holds for all x,y in X. 
with X < y, i.e. it then holds irrespective of the Birkhoff property of x and y. This uniform 



twist condition for instance holds for the Frenkel-Kontorova problem, see formula (1.51. 

To finish this section, let us for completeness include the following alternative existence 
proof for globally stationary Birkhoff solutions of arbitrary rotation vector. It was provided 
by Gole in [7] in dimension d — 1. The below is a more or less trivial generalization to higher 



dimensions, see also [TT]. As opposed to the results presented in Section 4.2 it also holds 
without Hypothesis C that requires that the Sj are coercive. The proof presented here is 
slightly shorter and more direct than the proof in [11] . 

Theorem 6.5. ^4/50 without the coercivity condition C, it holds that for every lu G M.'', there 
exists an X £ Blj with VW{x) — 0. 

Proof. Recall that the conditions A, D and E alone guarantee that the compact set Bu is 
forward invariant under the negative gradient flow. Condition B will be used below. 

Now, for B G Z'' a finite subset, recall the definition of the finite action Wb{x) := 
Ejgs 'S'i(^)- Then, for i G B^''\ it holds that diWaix) = diW{x), whereas if ||i - B|| := 
minjgs !1j ~ *11 > r, it is true that diWB{x) — 0. Thus, the time-derivative of Wb along 
solutions of ^ = —VW(x) equals 

J^WB{x) = ~Y.^^WBix)^^W{x) = ~ (9,W{x)f - ^ d,WB{x)diWix) . 

ieZ<i i6S('-) i^S('-),||i-S||<r 

We call ABix) ~ J2,eB^r-){d^W{x))^ . It is the square length of the gradient of the map 
y !-)■ Wb{x + y) from R^' ' to R evaluated at y = 0. With this definition, one checks 
that if B\ C B2, then AB-^^ix) < Ab2{x). Moreover, if B\ and B2 are disjoint, then 
AsiUB^ix) = Abi{x) + Ab2{x). 

The second sum in the expression for ^Wb consists of "boundary terms". We will call 
it aB{x) :— Yli<^Bir) diW B{x)diW {x) . Because is compact and diSj = di-jSo o tj^ for 
all j G Z'', there is a constant c > with the property that |i9iS'j| < c for all i,j G Z'^ and 
uniformly on B^. This in turn implies the estimate |aB(a;)| < c^(2r + l)'*|9i3|. 

Assume now that there is no globally stationary point in Bui. Then for every x G Bu, 
there is a finite subset Bx C Z'' such that Ab^{x) = 2ex > 0. Moreover, because VWb^ is 
continuous on R^ , it holds that x has an open neighborhood Ux C R^ on which Ab^ > Sx ■ 
By compactness we can find a finite collection xi, . . . , Xm G Bui such that Buj C yjYL\Uxi ■ 
Define B :— U]^iBxi and e :— min;{£a,, } > 0. Then every x £ Bui is in some Uxi and thus, 
Ab{x) > Ab^^ (x) > Exi > e > 0, that is Ab > e > uniformly on Bui- Moreover, translation 
invariance implies that for any G Z'*, also Ab+u > £ uniformly on Bui. 

For n G N, define the ball B{n) := {j G Z'' | \\i\\<n} CZ"^ and let iV G N be such that 
the B above is contained in B(N). Then Ab{n)(x) > Ab{x) > e > 0. 

Let m > 2 be an integer. By translation invariance and the fact that B{mN) contains 
at least m'^ translates of B{N) with disjoint r-interiors, it holds that ^s(mjv) > 'm'^AB(N) > 
m'^e. On the other hand, \amN{x)\ < c^(2r + l)''\dB{mN)\ = D'^m''-^ for some DJv > 0. 
Thus, -^WB{mN){x) < —m'^e + D^m''~^ and hence by choosing m large enough, we can 
arrange that ^WmN (x) < — 1 uniformly on Bui. 
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Since Bui is forward invariant under the negative gradient flow, this implies that WmN is 
not bounded from below on Buj- This contradicts the fact that Bui is compact and WmN is 
continuous. This proves that there must be a globally stationary point in Bui- □ 

7 Ghost circles 

In dimension d = 1, the concept of a ghost circle was introduced by Gole. We generalize 
this definition here to general dimensions. Note the similarity with Definition |4.13| of an 
Aubry-Mather set. 

Definition 7.1. [Ghost Circle] A ghost circle F C is a collection of configurations with 
the following properties 

• r is nonempty, closed and connected 

• r is strictly ordered, i.e. for every x,y£T,x<^y, x = y or x^y 

• r is invariant under shifts: if a; £ F, then for every {k, I) € Z'' x Z, also Tk,ix £ F 

• F is invariant under the positive and negative gradient flow: for all t £ R, 3't(F) — F 

An example of a ghost circle are the connected Aubry-Mather sets of Theorem |4.18[ 

The strict ordering and the shift-invariance of a ghost circle F imply that any configura- 
tion a; £ F is Birkhoff and hence has a rotation vector lj = ijj{x). The ordering of F moreover 
implies that this rotation vector is independent of the choice of x £ F, that is a; = tt'(F) and 
thus, T cBui. 

Let j £ Z''. Recall the definition of the projection to the j-th factor 

TTj : M.^'' -i> R, TTj{x) = Xj . 

Each TTj is continuous with respect to pointwise convergence. In fact, we can show that 
TTj |r : F — > R is a homeomorphism: 

Proposition 7.2. Let T be a ghost circle. Then, for every j £ Z'', the projection nj : R^ — >■ 
R induces a homeomorphism TTj\r : F — > R. 

Proof. Let F C R be a ghost circle. Clearly, 7rj|r is continuous. 

The strict ordering of F implies that nj\r is injective. Moreover, shift-invariance of F 
implies that if a; £ F, then so is tqjx — x + I for every I £ Z, whence the range of 7rj|r is 
unbounded. Since F is connected and 7rj|r is continuous, its range is both unbounded and 
connected, that is tTj It : F — >■ R is surjective. 

To prove that (tt^It)"^ : R — > F is continuous, it suffices to realize that nj : R^ — >■ R is 
an open map, i.e. that it sends open sets to open sets. This holds because the topology of 
pointwise convergence is generated by open sets 17 C R^ for which 7rfe((7) — V with C R 
an open subset, while ni{U) = R for all / 7^ k. For such U, it is clear that TVj(U) is open. □ 

Lemma |7.2| thus says that a ghost circle F is homeomorphic to R. It should be remarked 
though that, because F is invariant under the vertical shift tqj, and the gradient fiow 
'I't is equivariant with respect to ro.i, it makes sense to identify every element a; £ F with 
To,ix = a;+l £ F. The quotient F/Z ~ R/Z is a genuine topological circle. This identification 
is sometimes understood in this paper. 

The name ghost circle refers to the fact that F/Z may not consist of "physically relevant" 
configurations, i.e. globally stationary solutions. But, being a compact one-dimensional 
object consisting of orbits of a formal gradient flow, it has a good chance of containing such 
solutions. In fact, the following proposition serves as a first motivation to study ghost circles. 

Proposition 7.3. Every ghost circle F C R^ contains a globally stationary solution. 

Since F is a closed, fiow-invariant subset of some Bui, the proof of this proposition is identical 
to that of Theorem |6.5[ Moreover, we remark that when F contains at least one global 
minimizer, say x, then it automatically contains the entire Aubry-Mather set M{x). 

In the following two sections we first of all show that under generic conditions, ghost 
circles of rational rotation vectors exist and then we will prove a compactness result for 
ghost circles which will allow us to take limits and obtain ghost circles of irrational rotation 
vectors. 
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8 Morse approximations and periodic ghost circles 



In this section, we will prove two technical results. The first is that the local potentials Sj can 
be perturbed, in a way that will be made precise, so that the periodic action Wp^q : Xp,q — > R 
becomes a Morse function. 

The second result of this section says that whenever Wp^q : Xp,g — > R is a Morse function, 
then there exists a ghost circle F C ^p,q- 

Together with the results of Section |9] this will imply that any collection of local poten- 
tials admits a ghost circle of arbitrary rotation vector. 

8.1 Existence of Morse approximations 

Let uj € Q'* be a rational rotation vector and let (pi, qi), . . . , [pd, Qd) be a set of principal 
periods for oj. Recall that in Section[4]we defined the periodic action function Wp^q : Xp,, — >■ 

Rbyl¥p,,(x)=Ei6B,5.W- 

One says that Wp^q : Xp.^ — > R is a Morse function if at its critical points its Hessian is 

nondegenate. In other words, if \7Wp,q{x) = implies that D^Wp,q{x) is invertible, where 

D^Wp^q(x) is the symmetric matrix of second derivatives of Wp^q evaluated at x. By the 

implicit function theorem, every critical point a; of a Morse function is isolated. Moreover, 

each of these critical points can be assigned an index i{x) which equals the dimension of 

the unstable manifold of x, considered as an equilibrium point for the negative gradient flow 

^ = -^Wp,qix). 

We remark here that for arbitrary local potentials Sj, the periodic action Wp^q is not 
automatically a Morse function. A simple example of a non-Morse action function arises in 
the Frenkel-Kontorova model without local potential, for which 

Wp,q{^)= E ^ E (^«-^^)' ■ 

jeBp ||i-j||=i 

This action function satisfles Wp^q{x + t) — Wp^q{x) for all t £ R, so that its second derivative 
is everywhere degenerate. In fact, it has a one-parameter family of stationary points, and 
thus none of those is isolated. Nevertheless, in this subsection we will prove the following 
theorem: 

Theorem 8.1. Let Sj : R^ — >■ R be local potentials that satisfy conditions A-E. Let uj G Q'' 
and let {pi,qi), . . . , {pd,qd) he principal periods ui, that is Xi^ = Xp,,. Then there exists a 
sequence of local potentials S" with the following properties: 

1. The S" satisfy conditions A-E. 

2. The range of interaction of the S" is uniformly bounded in n. 

3. For every n, the periodic action Wp^q = X^jgB '•^ Morse function on Xp,,. 

4. The gradients converge uniformly: lim„^oo VS" = VSj uniformly in C^{M.^ ) 

5. The potentials converge uniformly on compacts: lim„_»oo — Sj uniformly on Bp^q. 

In dimension d = 1, this theorem was proved by Gole [8], [9] in the context of twist maps. His 
proof does not generalize to dimensions d > 1 or to general monotone variational problems 
in dimension d = 1, because it explicitly exploits the interpretation of Sj{x) — S(xj , Xj+i) 
as the generating function of a twist map of the annulus, see Appendix A. 

Our proof in higher dimensions is different, and it is based on Lemma |3.10| and ideas 
from equivariant Morse theory. We start by making Lemma [3.10| a bit more quantitative: 

Lemma 8.2. Letu) € Q''. Then the t -action of {U^ xX) / on^^ is properly discontinuous. 
More precisely, when (k, I) represents a nontrivial element of {JJ^ x 1)IL^ and x £ Xi^, then 

\Tk,ix - x\i := ^ \irk,ix-x)t\ > 1 . 

Proof. Let (p,q) be principal periods for uj, i.e. lj — —p^^q, and write n := | detp |. We 
notice that for an arbitrary G Z'* it holds that —n{k,uj) = (fc, |detp \p~'^q) £ Z and 
thus that {nk,—n{u),k)) € L^^. Hence, writing nl — —n{ui,k) -\-nl -\- n{uj,k), we see that 
Tk^iX = T„k,nix = To,„(i+(i^,fc))a;. Thus, |rfc_,a; - = \l + {(.d,k)\. 
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Now if {k, I) represents a nontrivial element of (Z'' x Z)/!^^, then n ■ \l + {lu, > 1, and 
hence we have that \tJ^iX — x\i > n. We claim that this implies that {t^jx — x\i > 1. This 

X = Tkfiijl iX — tI~i^x) and thus, by induction, that 



j+i _ j 

^k,l ^ '''k,l 



Wk'^i^^^'^k.M^ = \Tk,ix-x\x. Therefore, \tJ^iX~x\-i_ < [Tk^iX-rJ^^^ ^x\i + .. . + \Tk,ix~x\i = 
n\Tk,ix — x\i, which means that \Tk,ix — x\i > 1. □ 



8.2 



at hand, one can prove that the quotient Xtj/(Z x Z) is a smooth manifold 



With Lemma 

An arbitrary Z x Z-invariant function / : X„ — > R descends to this quotient and can hence be 
perturbed into a shift-invariant Morse function Instead of providing this rather standard 
construction from equivariant Morse theory, let us prove this latter fact directly here: 

Theorem 8.3. Let u) € Q'' and let p, q be principal periods for uj. When f : Xp,, -^M is an 
m > 2 times continuously differentiable shift-invariant function, then for every e > there 
exists a shift-mvariant Morse function f^ : Xp,g — >■ R with 

11/ - ,rilc"(B(iV)) < (J^^^ M every N >0 . 

Here, 

E{N) :— {x £ Xp,q I \xi — Xk\ > N for some i k with i,k £ Bp} . 
Proof. Let c := ^| ^^^^ and define the discrete collection of configurations 

Gp,q := {X : 7,''' ^ C ■ Z \ Tp-^q-X = X iOT j = I, . . . ,d } C Hp.q ■ 

We first of all remark that it is clear that Tk,iGp,q = Gp^q. For x £ Gp^q, let us now define 
the balls 

Br{x) ■— {y e Xp,5 I \y — x\i ■— ^ \xi - yi\ < r } . 

Then we have that Tk,iBr{x) = Brijk^ix), because \Tk.ix~Tk.iy\i = \x — y\\, that is the norm 
I ■ |i on Xp^q is shift-invariant. 

Moreover, if y G '^p,q; then there must be an element x G Gp,q with \xi — yi\ < |c for 
all i, that is for which |a; — y\i < |. In other words, Xp,, = U^gc Br{x) when r > |. 

implies that when r < |, then Br{x) n Brirkdx) = unless 
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On the other hand. Lemma 
{uj,k)+l = 0. 

This proves that for \ < t < \, the collection {Br{x)}x&Gp q forms a shift-invariant 
covering of Xp_q on which (Z"^ x Z)/I^ acts "properly discontinuously" . 

Finally, we let (j) : Xp,q — ^ [0, 1] be a C°° bump function with the properties that (j) = Q 
outside -Bi (0) and = 1 on _Bi (0). Let's say that | Ifl^l Ic^cxp ,) < E. 

After these preparations, we are ready to construct the perturbation of /. This is 
done by enumerating Gp,q = . . .} and defining it inductively. 

So let us assume that f^-i is r-invariant, satisfies the Morse property on the union 

Ui<i<n-i^i(2;') and fulfills the estimates ||/ - /^_i||c'"(B{iV)) < ""'^T+iv^)^ 

We now want a„ G R^^ to be a vector so that x i-> /^_i(a;) + {otn,x) is Morse on B i (cc"). 

4 

Such otn's are dense in R^^ by Sard's theorem, see for instance |10) . 
The function is now defined as the shift-invariant function 

fu{x)~f^--i{x)+ ^ ct){Tk,l{x ~ Xn)){a„,Tk,l{x - X^)) . 

(k,i)e{z<'xi,)/i^ 

Because B i{Tk,iXn) does not intersect B i{TK,LXn) unless {k,l) — {K,L) mod I^j, we have 
that at every x G Xp^,, the above sum consists of only one term. Moreover, is Morse on 
Bi{x^) by construction. 

In fact, by choosing a„ small enough, one can make sure that is Morse on the entire 
union Ui<i<n Bi{x'^). This is true because the collection of Morse functions is open in the 

space of differentiable functions C™ (^Ui<i<n-i ^'^^ Tn>2, see [10) . 

By choosing Q!„ even smaller if necessary, we can also arrange that fn-\ — fn has a 
C""(i5(A''))-norm less than j^^p^^fp . This implies that 

(l-2-("-^))£ 2-"£ ^ (l-2-")£ 
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The required is the hmit /' := hm„_>cx) fn- Not only does this limit satisfy the required 
estimates, but it also stabilizes pointwise, which shows that it is Morse. □ 

We can now complete the proof of Theorem |8.1| 



Proof of Theorem \8.1\ We start by perturbing the Sj so that they satisfy a strict monotoni- 
city criterion. This will then allow us to perturb the potentials once more without risking to 
destroy monotonicity condition D. Recall that Bp ~ p([0, l]**) nZ'' is a fundamental domain 
of p. It has cardinality | detp |. Our first perturbation step is now made by defining 

Sf{x) := Sj{x) H ^ {xk — Xi) arctan(xfc — Xi) . 

The strict monotonicity of the S" follows because (xarctana;) = (-^^^a-^a is strictly 
positive. Hence, 

di^kS"(x) < -J- — — < for all i.k £ j + Bp with i k . 

n (1 + (xi - Xfe)^)^ 

the periodic action Wp,q : Xp,, — ^ R defined by Wp^g{x) := J2jeBp 
can now be perturbed into a r-invariant Morse function Wp^q : Xp,q — > R of the form 

W^,q{x) = W^,q{x) + F-{x) . 

The perturbation may be chosen so that it satisfies ii-F"|jc2(Xp ,) < F"{'''k,ix) = F"(x) 
for all X £ Xp,q and all k,l and \di,kF"{x)\ < ^ (i+(^.i^^)a)2 ■ 

For a; : Z"* — ^ R, let us denote by x\'^'^g the p-periodic extension of x\j+Bp defined by 
(^1 j+s ) . = Xk, where k £ j -\- Bp is the unique element of j + Bp equal to i modulo p(Z'') . 
Then we can define, for each j G the new local potential 

Snx):-~Slix) + -^F-{x^-,^) ■ 



By Theorem 
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We will now prove that these satisfy all requirements of Theorem 8 

In fact, condition A and requirement 2. hold true because the range of interaction of 
both the sum ^ "Y^^ kej+B (^^ ~ -^0 arctan(2;fc — Xi) and the perturbation F" (^x\'^^g ^ do 
not exceed the bounded radius of Bp. 

Condition B holds by definition. Condition C holds because x i-^ x arctan x is non- 
negative and \F„{x)\ is uniformly bounded. Condition D holds true because di,kS^{x) < 
di^kSjix), as is easy to check. 

Requirement 3. holds because = Yj^B ^ji^) ^ Wp,q{x) + F"{x) is a Morse 

function by construction. 

Requirement 4. and condition E are true because both | ^ a; arctan x] — | j_^^'2 +arctanx| < 

3 and | a: arctan x| = | i^i_^_^^2p I ^ 2 are uniformly bounded and ||-F"||02(Xp < 2n' ^'^ that 
lIVS" — VS'jll^ijjj-dj < ^ for some constant C depending on the dimension d, the period- 
icity p and the range of interaction r. 

Similarly, requirement 5. holds true because kej+B ~ •'''=) arctan(a;i — Xk) is uni- 
formly bounded on Bp^q and |F„(a;)j < uniformly on Xp,,. □ 



8.2 Existence of periodic ghost circles for Morse actions 

We will now show that when the local potentials Sj satisfy conditions A-E and are chosen so 
that Wp^q : Xp,q — ^ R is a Morse function, then they admit a periodic ghost circle F C Xp,q. 
More precisely, we will prove the following: 

Theorem 8.4. Let to G Q'' and let (pi, gi), . . . , [pd, qd) he principal periods for uj. Assume 
moreover that the local potentials Sj are chosen so that Wp^q : Xp,, — >■ R is a Morse function. 
Then there exists a ghost circle F C Xp,, for the Sj . This ghost circle includes all 
the global minimizers of Wp^q. It consists of stationary points of index and index 1 and 
heteroclinic orbits of the negative gradient flow. 
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The construction of this ghost circle is essentially the same as the construction in dimension 
d = 1 provided by Gole We nevertheless decided to provide the proofs. 
To prove Theorem |8.4[ we need two lemmas and the following definition: 

Definition 8.5. We say that x <gi y are consecutive index-0 stationary configurations if 
there is no index-0 stationary configuration z with x <^ z <^ y. 

It turns out that when Wp^q is a Morse function, then between consecutive index-0 stationary 
configurations we can find another critical point: 

Lemma 8.6 (Mountain pass theorem). Assume that Wp,q : Xp., — >■ R is Morse and let 
X y be two consecutive index-0 stationary configurations of —\/Wp,q. Then there is an 
index-1 stationary configuration z in between x and y. 

Proof. We use a simple variant of the mountain pass theorem, see for instance [H], Section 
8.5.1. For this purpose, we let C be the collection of curves from x to y lying in the order 
interval K := [x,y], that is 

C := {7 : [0, 1] — )■ if I 7(0) — X, 7(1) — y and 7 is continuous} . 

We now claim that there is a critical point z £ K for which Wp^q(z) — c, where 

c :— inf max Wp ailit)) ■ 
7GC0<t<l '^'^^ 

To prove our claim, let us define, for 5 £ R, the sub-levelsets 

--{xeK \ Wp,q{x) < S} . 

These are invariant under the forward flow of ^ = -'\/Wp^q{x). This is true because K 
is invariant and because Wp^q is a Lyapunov function for the gradient flow. 

Suppose now that there is no critical point x <^ z <^ y with Wp^q{z) — c. We will show 
that this leads to a contradiction. 

We first of all remark that, by the Morse lemma and the fact that x and y have index 0, 
it holds that c > max{Wp,g(a;), Wp^q{y)}. Thus, because there are only finitely many critical 
points in K, and none of these except x and y lie in dK, there exists an e > so that the 
set K''^^\K'^~^ does not contain any critical points. 

This in turn implies, by compactness, that there is a cr > so that ||VWj,,q||^ > a 
on K'^'^^\K''~^^^ . Hence, a solution curve t 1— >■ x{t) of the negative gradient flow satisfies 
£Wp,q{x{t)) = ~\\\/Wp.qix{t))\\'^ < -a SO long as x{t) G K^+^XK"-"^^ . In particular, there 
is a T > for which *t(A"=+=) C K'"''^^. 

At the same time, by definition of c, there exists a 7 G C with 7([0, 1]) C K'^'^^ . The 
curve o 7 G C then lies entirely in K''~°^^ . This contradicts the definition of c and hence 
there must be critical points x ^ zi, . . . , Zm ^ y with Wp,q{zi) — c. 

It remains to show that at least one of the Zi has index one. In fact, the argument is a 
bit subtle. We start by observing the following: 

1. \i X <^ Zi <^ y is an index-0 critical point with Wp,q{zi) = c, then there are ai,/3i > 
so that whenever 7 G C intersects Bc^ (z), then maxtg[o,i] W'p,q(7(i)) > c + 

2. If 2; <C <C 2; is an index-> 2 critical point with Wp^qizC) = c, then there is an > 
so that whenever 7 G C intersects Baiizi), then 7 is homotopic to a curve 7 G C 
with the property that 7 does not intersect Ba^^Zi), while maxtg[o,i] Wj,,q(7(i)) < 
maxts[o.i] Wp,g(7(t)). 

These statements are easy to prove in local Morse coordinates near the critical point Zi. 
At the same time, by compactness, we have that there exist e,ai,(T2 > so that a\ < 
||VVyp.,|P < a2 on uB;3i/2(«i) U . . . U B^„/2(z™)) . Using that 1 1 1 1 = 

||VWp.g(x(f))|| and -^Wp.q{x(t)) = —\\\/Wp^q(x{t))\\^ for solutions of the gradient fiow, one 
can prove quite easily that this implies: 

3. If t i-> x{t) solves ^ = --VWp,q{x) and Wp,q{x{0)) < c 4- e and a:;(0) ^ Bi3^{zi)LI 
. . . U Bi3^{zm), then for < t < T := min{/3i, . . . ,/3m}/2^/(T^ one has that x{t) ^ 
U™ iB^,/2(zi) and hence Wp,q{x{T)) < max{c - e/2, Wp,q{x{0)) - Tai}. 
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We now use these facts as follows: Let us assume that none of the Zi has index 1 and let 
7n £ C be a sequence of curves with c < maxtg[o,i] Wp^qi^nit)) < c + ^. By property 1. 
we know that for large enough n, the curve 7„ does not intersects Bon{zi) for any of the 
index-0 points Zi. At the same time, by property 2. we may assume that none of the 7n 
intersects the Ba^i^Zi) for any of the index-> 2 points Zi. Property 3. then implies that for 
large enough n we have that (v&t o 7n)([0, 1]) C K'^~^ for some 5 > 0. This contradicts the 
definition of c. □ 

The next step is to show that the unstable manifold of the index- 1 critical point z of Lemma 
defines ordered heteroclinic connections to its neighboring index-0 critical points x 



.6 



and y. This result, in more generality, can also be found in [T], see Theorem 1 of Chapter 4. 

Lemma 8.7. Let Wp^q : Xp,, — >■ R 6e Morse, let x <^y be two consecutive index-0 stationary 
configurations ofWp^q and let z be an index-1 stationary configuration with x <^ z <^ y. Then 
the unstable manifold of z forms strictly ordered heteroclinic connections from z to x and y. 

Proof. We consider the linearization of the negative gradient vector field at z, given by 
the matrix —D^Wp^q{z). The twist condition D and the bound on the second derivatives 
E together guarantee that there exists a constant M > so that the symmetric matrix 
H := -~D^Wp,q{z) + Mid is nonnegative and strictly positive on its diagonal and its two 
off-diagonals. 

By the theorem of Perron- Frobenius, H then has to have a unique simple largest eigen- 
value Amax + M € R+ and the corresponding eigenvector Cmax can be chosen strictly positive. 
Because z is an index-1 point, Amax is then the unique positive eigenvalue of —D'^Wp^q{z) 
and emax is its strictly positive eigenvector. 

The unstable manifold W"(z) of z is one-dimensional and at z it is tangent to emax. In 
fact, it consists of z and two orbits of the negative gradient flow 



a" 



~^{t) ^ z± e^"'°'''* emax + o(e^"'"'*) for t — ^ — oo 



In particular we see that close to z, the unstable manifold is strictly ordered, because emax is 
strictly positive. Theorem |6.2| then implies that the entire W"(z) is strictly ordered. Thus, 
we see that there must be two critical points z~ := limi_>oo 0-{t) and z^ :— limt^oo a+(t). 
We claim that z~ = x and z^ — y. 

To prove this, we will show that z~ and z'^ are index-0 critical points. Our claim then 
follows because x < z~ ^ z <^ z^ < y hy monotonicity and because x and y are consecutive 
index-0 points. 

So let us consider the linearization matrix -'D'^Wp^q{z~). It also has a unique maximal 
eigenvalue Amax a^nd positive eigenvector e'ax- We know that limt_»oo ce~it) = z~ and that 
a~(R) is strictly ordered. At the same time, because —D^Wp,q{z~) is symmetric, emax is 
perpendicular to all other eigenvectors of —D^Wp^q{z~), which implies that none of these 
other eigenvectors lies in the positive or the negative quadrant. This means that a~ {t) has 
to approach z~ tangent to emaxj that is 

ce~ {t) = z~ + e'^'"'"'*Cmax + o{e'^'"^^*) for t ^ oo . 

In particular, Amax < 0. But Amax is the maximal eigenvalue of —D^Wp.q{z~). This means 
that all eigenvalues of D'^Wp^q{z~) are positive, i.e. that z~ is an index-0 point. 

A similar argument for finishes the proof. □ 

We conclude with a definition and then give the proof of Theorem |8.4[ 
Definition 8.8. A nonempty and strictly ordered collection of configurations 



Co = {..., 2;_i < a;o < a;i <...} C 

x-Q skeleton for the Mo 
it consists of index-0 critical points of Wp^q 



is called a maximal index-0 skeleton for the Morse function Wp,q if: 



• it is shift-invariant: for all x £ Co and (fc, I) £ x Z, it holds that t^^ix £ Co 

• it is maximal: if y ^ Co is an index-0 point, then there is no i G Z with Xi <^y <^ Xi+i. 
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Proof of Theorem \8.4\ We remark that a maximal index-0 skeleton in general is not unique, 
but it is not hard to see that a maximal index-0 skeleton exists if Wp^q is Morse. 

Indeed, one can construct one by starting with the strictly ordered, shift-invariant col- 
lection Co = {. . . , X-i,xo, xi, . . .} of all the global minimizers of Wp^q. We note that Co is 
discrete because Wp^q is a Morse function. 

If there exists an index-0 point x ^ Cq with the property that Xi x <^ Xi for some 
i, then one augments Co by the r-orbit of this x, thus obtaining the strictly ordered, shift- 
invariant and discrete collection 

Co' := Co° U {Tk,ix I (fc, l)eZ''xZ} . 

One keeps on adding r-orbits of index-0 points this way. The Morse property of Wp^q guar- 
antees that the number of index-0 points between xq and + 1 is finite, which implies that 
this process stops after finitely many steps. 

The maximality of an index-0 maximal skeleton Co = {. . . , x_i <^ xo <^ xi <^ . . .} just 
means that the pairs Xi,Xi+i are consecutive index-0 points. The Mountain Pass Lemma 
guarantees that between these consecutive elements, there is an index-1 critical point Zi, 
while Lemma |8.7| says that the unstable manifold of this Zi defines strictly ordered hetero- 
clinic connections from Zi to Xi and Xi+i. 

If we choose the Zi in such a way that Co := {. . . ,a;_i <C z-i <ti xq <^ zq <^ Xi <^ 
z\ <^ . . .} is shift-invariant, then the union of Co and these heteroclinic connections is a 
ghost circle F. The construction above shows that F may be assumed to contain all global 
minimizers of Wp,,. 

It only remains to show that this F is C'. This is clear except at the critical points. 
But in the proof of the Lemma [8. 7[ we have seen that at the critical points, the heteroclinic 
connections are tangent to the dominant eigenvector. This eigenvector is simple and hence, 
r is C' also at critical points. □ 



9 Convergence of ghost circles 

Section |8] was devoted to the construction of periodic ghost circles for action functions that 
satisfy the Morse property. In this section we will prove the existence of periodic ghost 
circles for arbitrary action functions. In turn, this will then imply the existence of ghost 
circles with irrational rotation vectors. These results follow from a compactness theorem for 
ghost circles that we will prove below. Before we can formulate it, let us specify what it 
means for a sequence of ghost circles to converge: 

Definition 9.1. [Convergence of ghost circles] We say that a sequence of ghost circles r„ 
converges to a ghost circle Foo , if for every ^ £ K, the sequence of configurations a;" (^) G F„ 
defined by 7ro(a;"(^)) — ^ converges pointwise to the configuration x°°{^) G Foo defined by 

Thus, if T^n ^ Foo as n )■ oo then Foo consists of pointwise limits of elements of the F^. 

Before stating the most important results of this section, let us make a few simple ob- 
servations concerning convergence of ghost circles. First of all, one can observe that if 
lim„_>oo i^n ~ ij^oo and if a sequence of ghost circles F,i C Bui„ converges to Foo, then it 
must be true that Foo C Bui^ ■ This follows from the continuity of the rotation vector as a 
function on B and the fact that the rotation vector of a ghost circle is defined as the rotation 
vector of any of its elements. 

The second remark is that if the F„ C Bu) are periodic ghost circles with the same ratio- 
nal rotation vector, and F„ — > Foo, then Foo C Bu) is periodic as well. This follows because 
Buj is a closed subset of . 

Our compactness result now is the following: 

Theorem 9.2. Let ujn €z K C R"* be a sequence of rotation vectors contained in a compact 
set K and converging to ujao € K C M.'' and let 5" : M. be a sequence of local potentials 

such that VS" converge to VSj^ uniformly in C"'^(R^ ). Finally, let F„ be a sequence of ghost 
circles for the S" of rotation vector cj„. Then there exists a ghost circle Foo for the 5°° of 
rotation vector ujaa and a subsequence {nj}|5ii such that limj^oo F„j = Foo. 

// moreover lim„_>oo 5*" = Sf^ uniformly in C^\Bk) and if every Fn contains a global 
minimizer, then also Foo contains a global minimizer. 
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Before proving this compactness result, let us formulate its two main implications: 

Theorem 9.3. Let uj G Q** be arbitrary and let the local potentials Sj be given. Then there 
exists a periodic ghost circle Fi^ C Bu for the Sj . This Tuj may be chosen so that it contains 
a global minimizer. 

Proof. Given lj G Q'' and any local potentials Sj, choose principal periods {p,q) for cj. By 
Theorem 



8.1 



we can choose a sequence of local potentials S" such that lim„^oo — ^■S'j 
uniformly in C^(R^ ) and lim„_>oo S" = Sj uniformly in C'^{Bk), while at the same time 
Wp^g : Xp,q — > R is a Morse function. Then, by Theorem 8.4 there is a ghost circle T n C 



Bp^q = Bui for the local potentials 5*" that contains a minimizer of Wp^g. By Theorem 9.2 
a subsequence of the r„ converges to a ghost circle F,^ C Bp^g = B^i for the local potentials 
Sj. By the second conclusion of Theorem 19.21 F^j contains a global minimizer. □ 



Theorem 9.4. Let lj G IR''\Q'' and let the local potentials Sj be given. Then there exists 
a ghost circle F^j C Bui for the Sj. This Fui may be chosen so that it contains the entire 
Aubry-Mather set of rotation vector u. 

Proof. Given uj £ and local potentials 5*-,, choose a sequence cj^i G Q'' such that 



lim„_>txj a;^ = cj. By Theorem 9.3 there is a periodic ghost circle r„ C Bui„ for the local 
potentials Sj that contains at least one global minimizer. By Theorem |9.2[ a subsequence 
of the r„ converges to a ghost circle Fui C Bui- 

The requirement for the second conclusion of Theorem |9.2| is trivially valid, so that Fui 
contains a global minimizer, say x. Being closed and shift-invariant, this implies that Fui 
contains the entire Aubry-Mather set Mix). □ 



Before proving Theorem |9.2[ we remark that if r„ is an arbitrary sequence of ghost circles 
for the local potentials 5" and with rotation vectors aj„ in a compact set K, then for every 
^ G R the sequence of configurations a;"(f) G r„ has a subsequence that converges pointwise. 

This just follows from the compactness of Bk H G R* j 7ro(a;) = The problem is to 
show that this subsequence can be chosen independent of ^ and that the collection of limit 
configurations {lim„_»oo a^"(C) I C G 1^} forms a ghost circle for the Sj — lim„^oo S". 



We will now make some preparations for the proof of Theorem |9.2| To start with, we define 
for a given ghost circle F, the map 



F by - *_io(^o|r)"\ that is: T^{^)~ {t"^ iO) = K o o K|r)"') (0 
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Here, *_i : X -i> X denotes the time— 1 flow of |^ = -\/W{x). By Theorem 
homeomorphism, being the composition of two homeomorphisms. Moreover, it is "pointwise 
Lipschitz continuous" : 

Lemma 9.5. Let K C R'' be a compact set and F C Bk = ^uii=kBui a ghost circle with 
rotation vector uj G K for the local potentials Sj satisfying conditions A-E. Then, for every 
k G Z"^, there is a constant A||fe|| > 0, depending only on K and jjfcjj such that 

\Tl{i)-Tl{y)\<k\\k\\\i-i^\. 

Proof. Let ^, G R and denote X — (7ro|r)^^(0 ^-iid Y = {-roIy)'^ [v). Assume that <u, 
whence X < y. Denote by *t the time-t flow of -VW^, with W ~ Y.jez'^ Sj. Since F 
is forward and backward invariant under 'I', both ^'_i(X) and lie in F and satisfy 

^_i(X) <C Now we apply the parabolic Ifarnack inequality of Theorem 6.4 to 

t = 1, i — 0, X = and y — 'i/-i{Y), to find that there is an L > depending only 

on K and such that 

Tfc M - T^iO = - <l(Yo- Xo) ^li^^-O ■ 

A similar argument in the case that ^ > finishes the proof. □ 

We remark here that we see no reason why the maps TTk o (ttoIt)^^ should be uniformly 
Lipschitz continuous. This is why we study the maps o vp-i o (vrolr)^^ instead. 
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Definition 9.6. We say that a sequence of maps T" : M — converges point-wise uni- 
formly to a map r°° : R — )■ as n — >■ oo if for every k € Z"^ the sequence of maps 
Tfc" := TTfc o T" : R "> R converges uniformly to Tfe°° — tt^ o r°° as n oo . 

Corollary 9.7. Let K G M.'^ be a compact set. Assume that for every n € N, we are given 
a rotation vector ujn G K, local potentials S" satisfying conditions A-E and ghost circles 
Tn C Bui„ for the local potentials S" . 

Then there is a subsequence {nj}j^E C N with the property that the maps T^^j : R — >• R 
converge pointwise uniformly on R, say T^^j — >. T°° as j — >■ oo. Each limit map := 
TTfe o T°° : R — ^ R, is non- decreasing, surjective and Lipschitz continuous. 



Proof. Fix a fc £ Z"*. By Lemma 9.5 the maps T^" :— vrfcoT'"" are uniformly-in-n Lipschitz 
continuous with Lipschitz constant A[n.||. Moreover, by the definition of a ghost circle they 
are 1-periodic and increasing. To see th at th ey are uniformly bounded on compacts, we 
then just have to note that Proposition 7.3 implies that Tq{0) £ [—1,1], while Lemma 



3.5 



and the fact that o (ttq) is a Birkhoff sequence then imply that T[{0) 

\K\\ ■ ||A;||,2 + \\K\\ ■ \\k\\], where ||iC|| := max^gx 
Thus, the theorem of Arzela-Ascoli guarantees that there exists a uniformly convergent 
subsequence Tj. — >■ T^^ for j — > oo. Clearly, Tj^ is nondecreasing and Lipschitz contiimous 
with Lipschitz constant A|H;|[. 

Let j H-)- fcj , N — >■ Z'' be a denumeration of Z'*. Then the diagonal sequence {njjjgN C N 
defined by Uj :— nj^kj has the property that T^"i — > T°° pointwise uniformly as j — >■ oo. □ 

Theorem 9.8 (Convergence of ghost circles). Let T„ be a sequence of ghost circles for the 
local potentials S" . Assume that there are local potential functions such that VSf 

VS°° uniformly m C^(R* ) and that the maps T^" converge pointwise uniformly. Then 
there is a ghost circle Too .for the local potentials S]^ such that r„ — > Too as n oo. 

Moreover, when r„ contains a global minimizer x„ and lim„_j.oo 5*" = uniformly in 
C^{Bk), then Too contains a global minimizer as well. 

Theorem |9.2| now follows directly from Corollary |9.7| and Theorem |9.8| 

Before we prove Theorem |9.8| let us recall that if Too = lim„_ioo r„ exists, then it must 
be equal to 

Loo := ■= lim ^"iO pointwise | C G K} • 

rt—^oo 

At this point, it is of course not clear whether the limit lim„^oo exists for every ^ £ R. 

To see that it does under the conditions of Theorem [9^ we note that = *i 

so that 

lim = lim *i (r"(C)) exists and is equal to ^^(^""(O)- 



This is true because on the one hand, according to Corollary |5.9[ v!/" — )■ uniformly in 
the topology of pointwise convergence, while on the other hand it holds that for every ^ £ R, 
the sequence of configurations T'""(^) £ B converges pointwise to the configuration T'°°(^) 
as n — > oo, because T^" — > T°° pointwise uniformly. Thus we find that under the conditions 
of Theorem |9.8[ Too '■= lini„^oo r„ C B:^^ is well defined and moreover that, if Too is a 
ghost circle, then T^-^ = T°° . 

We will now show that Too is in fact a ghost circle for the local potentials S^: 

Proof of Theorem \9.8\ We first check that Too has the properties required for a ghost circle: 

1. Closedness: Let a;°°(^m) £ T'oo be a sequence of configurations that converges point- 
wise. This implies that the converge, say to ^. We now want to show that 
lim„i->oo a;°°(Cm) = £ Too pointwise. This follows because limm^oo {^m) = 
lim„,^oo*?°(T°°(Cm)) = *r(lim„^oor°°(U)) = W{T°°iO) = AH these 
limits are pointwise. We have used that ^PJ" is continuous for pointwise convergence 
and that lim„_+oo r°°(Cm) = T°°(^) pointwise. 

2. Connectedness: We note that Too = $i°(T°°(R)), so it is the image under a conti- 
nuous map of a connected set, hence connected. 

3. Strict ordering: Suppose < v. Recall that T°° is nondecreasing, so T°°{£^) < 
T°°{i'). We remark that r°°(^) cannot equal T°°{v), because this would imply that 
^ = 7ro(2;°°(0) = 7ro(*r(r°°(C))) = MW{T°°{u)) = ^o(a;°°H) = v. Thus r°°(C) < 

r°°(i/). 
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The strict monotonicity of the negative gradient flow, then irnphes that x°°{^) = 

*j°(T°°(c)) < = x°°{u). 

4. Shift-invariance: Let x°° G Too, that is 2;°° = hm„^oo x" with x" e r„ and no{x") = 
Tvo{x°°). Let k £ and Z G Z be given. We want to show that Tk,ix°° € Too, that is we 
want to show that Tk^ix°° — hm„_j.oo y" with y" £ r„ such that 7ro(2/") = 7ro(Tfe_;a::°°). 
We prove this by writing 

lim (rfc,iX - J/ ) = lim (jk^ix - Tk,iX ) + lim (rfc^jx - t/ ) 

n— foe n — J-oo n — foe 

and showing that both hmits on the right hand side vanish. 

The first limit is zero because, by Lemma [5.3[ Tk.i is continuous in the topology of point- 
wise convergence. Thus we have that hm„_).oo {Tk.ix°° — Tk,ix") = lim„^oo Tk,o{x^ — 
a;") = 0. 

For the second limit, we realize that Tk,ix" £ r„ because r„ is shift-invariant and we 
observe that Tvo{rk,ix") = x^ + I. Because T^" — > r°° pointwise uniformly, we more- 
over know that lim„_i.oo T^"{xl + 1) = T°°(a;^ -I- I) pointwise. Thus, by the uniform 
convergence of the to VPf", 

lim {Tk,ix")^ lim *'i'(T"(4'-fO) = *r(r°°(a;r+0)= lim {T" [x^ +1)) = limy". 

ri— foo n— foo n— foo n—^oc 

5. Flow-invariance: This is proved in a similar way as shift-invariance. So, let x°° £ Too, 
that is x°° = lim„^oo x"^ with x" G r„ and tto{x") = 7Vo{x°°). Let t G R be given. We 
want to show that 5'^x°° G Too, that is we want to show that ^PJ^a;"" — lim„_>oo y" 
with y" G r„ such that 7ro(y") = TTo{'^rx°°). We prove this by writing 

lim - y") = lim - *"a;") + lim {^['x" - y") 

n— foo n — foe n — foe 

and showing that both limits on the right hand side vanish. 

The first limit is zero because, by Theorem |5.4[ converges to uniformly. Thus 
we have that hm„_>oo *ra:" = "i^x"". 

For the second limit, we realize that ^'"x" G P" because r„ is flow-invariant and we 
observe that lim„^oo 7ro(^'"2:") — 'K(){'i''^ {x°°)) because tto is continuous for point- 
wise convergence. Because T'"" — T°° pointwise uniformly, we therefore know that 
lim„^oer^''(7ro(*r(a:"))) = T"^ iM'^Tx"°)) pointwise. Thus, by the uniform con- 
vergence of the to "if, 

lim ^^2;"= lim *l'(r^"(7ro(*ra;")))= lim ^rC?''" (7ro(*ra;°°))) = lim y". 

n— foo n — foe n — foe n— foe 

We finish the proof of Theorem |9.8| by proving that when each r„ contains a minimizer and 
lim„_foe = SJ" uniformly in C^'(Bk), then also Too contains a minimizer: 

Minimizing property: Suppose that every ghost circle r„ contains a minimizer x" = 
x" This means that for every finite subset B C Z'' and every j/ : Z"* — >■ R with finite 
support in i?''^' it holds that 

WS{x" +y)-WS{x'')>0 , where WSix) ■- ^ S"(x) . (9.14) 

ies 

By compactness of Bk/Z, a subsequence of the x"'{^„) converges pointwise, say to x°° = 
limj^oe x "^ {(,„, )■ Moreover, lim„_»oe = Wb uniformly in C°{Bk)- Taking the limit of 
equation ( |9.14[ ) as Uj -> 00 then shows that W^{x°° + y) — Wg{x°°) > 0. In other words, 
x°° is a global minimizer. 

It remains to prove that x°° G Foe. This holds because x°° = limj^oo 2:"^ (^„^ ) = 
lim,_+oe *?(r"(C„^.)) = *r(r°°(Coe)) = a;°°(^oe), where ^oe - lim,_^oeCn,. □ 



10 Gap solutions 

In this final section we examine the situation that an Aubry-Mather set M{x) C Bu, has 
a gap, that is when there are elements y~ ,y^ G M{x) with y~ <^ y'^ such that [y~,y'^] 
does not contain any elements of A4{x) other than y~ and y'^ . This situation occurs when 
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oj G or when cj G R^XQ" and M{x) is a Ca ntor s et 



10.7 



below, which states that either [y , y'^ 



The main result of this section is Theorem 
admits a foliation by global minimizers, or there exists at least one stationary configuration 
z G [y~ ,y'^] that is not a global minimizer. This result is more precise than the result of 
[3], that says that a gap must contain at least one stationary solution. Moreover, the proof 
below is more geometric, as it makes use of ghost circles. 

We start with the following theorem, which is a refinement of a result by Moser [15]. It 
says that when a gap admits a foliation by stationary points, then all of them are minimizing. 
Recall that every Aubry-Mather set is contained in a ghost circle. 

Theorem 10.1. Let [y~ ,y'^] be a gap in the Aubry-Mather set M{x) and let V be a ghost 

circle so that M{x) C F. If F^^ F n [y^,y^] consists of stationary configurations 

only, then all of them are global minimizers. 

Proof. Assume that F'" '^^^ consists of stationary points only but that w G F'*' is not 
a global minimizer. Then there is a finite subset B <Zl/' and a ^ G with supp(z) C -B'""' 
such that Wsiyj + z') < Wb{w). Because the function z i— >■ Wb{w + z) is coercive, it attains 
its minimum, let's say at a Z with support in B^^\ By assumption Z 0. Let's say there 
IS an I G B'''' for which Z, > 0. In the case that Zi < the proof is similar. We now claim 
that Z can be chosen so that w + Z < y^ . 

To prove this claim, we remark that when m := min{m + Z,y'^} and Af :— maxj io + 
Z,y+}, then Wb{w + Z) + Wfl(2/+) > WBim) + Wb{M), as in the proof of Lemma liT 



Because both {w + Z) — m and y^ — M are supported in i?''^' and both w + Z and y 
minimize Wb with respect to variations supported in B^-^\ it must therefore hold that 
WBiw + Z) = Wsim) = Wb{w + min{Z,y+ ~w}). 

The next step is to define y :— mi{y (zr\y^w + Z}. Because w,y,y'^€r,w + Z<y^ 
and Zi > 0, it now holds that w <^ y < y'^ . At the same time, because F is connected, y 
touches w + Z. That is: there is an i G -B'*^' so that yi — Wi + Zi, while Zk + Wk ~ Zk < yk 
for all k ^ i?''"'. We claim that this is impossible. 

To prove this, choose such an i G -B'-'^' at which yi = Wi + Zi and a k G Z'' for which 
|i — fc|| = 1. Then, because y is a global stationary point and w + Z is stationary for Wb 
with respect to variations in B''^\ it must be true that 



= d,WB{y)-d,WB{w + Z)= [1 9,,iS,{ty + (1 - t){w + Z)) dt y{yi~w,^ Z,) 

jeB,iei'^ 







: / di.kSi 
Jo 



[ty + (1 - t){w + Z)) dt ■ (yk - Wk - Zk) 



Here, the inequality holds because yi — Wi — Zi — and di,iSj < when i ^ I and yi — wi — Zi > 
for all I. The twist condition that di,kSi < then guarantees that yk = Wk + Zk- By 
induction, one then finds that there is a ^ S''"' for which yt = Wk + Zk = Wk. This is a 
contradiction. □ 

We will now show that when F''' does not consist of only stationary points, i.e. mini- 
mizers, then it contains at least one non-minimizing stationary point. We do this by finding 
a stationary point of a "renormalized action" function : [y~,y^] — ^ [0,cx3). In or- 

der to define W^y- y+-^, we need the following well-known technical result that states, when 
applied to M — A4{x), that the gaps of an Aubry-Mather set are uniformly summable: 

Theorem 10.2. Let Ml be any strictly ordered, shift-invariant collection of configurations 
of rotation vector a; G R''. Let x,y £ M be so that x <^ y and assume that there exists no 
z e M with X <s: z ^y. Denote H^i ~ {i e Z"^ \ {oj,i) G Z}. Then 

Proof. We start by remarking that our assumptions on Ml imply that x and y are Birkhoff. 
Now, let i and j be representatives of different equivalent classes of 11^ jPL^ and let Tn,n £ Z 



be arbitrary. Then, (uj, j — i)-{-m—n 7^ and hence, by Proposition 3.8 either Tj-i^m-ny y 
or Tj-i^rn-nX 3> X. lu the first case, actually Tj-i^rn-ny < x because there is no element of 
M between x and y. Evaluating the latter inequality at i, we then obtain j/j + m < -f n. 
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In the second case, one finds that rj-i^rn-nX > y, whence xj + m > yi + n. In both cases we 
find that {xj + m, yj + m) n {xi + n,yi + n) — III. 

For 5 G R, denote by [^] := max{n G Z | n < ^} and define x — mii{xi — [xi]) G [0, 1). 
Then, clearly Xi — [xi] > x. We claim that j/i — [xi] < x + 1. To prove this, note that Xi — [xi] < 
Xj — [xj] + 1 and hence, by the Birkhofi' property, x <^ Tj-iX + [xi] — [xj] + 1. The assumption 
that X and y are consecutive elements of A4 then implies that y < Tj-iX + [xi] — [xj] + 1, 
that is, yi — [xi] < Xj — [xj] + 1. Hence, yi — [xi] < x + 1. 

This yields, denoting by \A\ the Lebesgue measure of a set A C R: 



i£Z'i/H^ iel.''/H^ 



[J {X^ - [Xi],yi - [Xi] 
ieZ<i/H^ 



< \{x,x + l)\ = 1 



□ 



For rationally independent rotation vectors, for which H^j = {0}, Theorem 10.2 was stated 
for the first time by Moser |14| . 

For a gap [y^ ,y^], with y~ ,y~^ £ Bui, let us define 



[y ,y+]:={yG[y ,y ]\ rk.iy = y iHi^,k) + i = o}. 



It is not hard to see that C Buj- Namely, when y^ < y < y and {uj,k) -\-l > Q, 

then Tfcjy > Tk^iy~ > y^ > y, where the second inequality holds because [y~,y^] is a gap. 
Similarly, Tk^iy < y when {u, k) + I < 0. Hence, y is Birkhoff once r^^iy = y for all k, I with 
{uj, k) + I — 0. We are now ready to define the renormalizcd action function; 



Definition 10.3. When [y ,y] is a gap, we define W^y- : — >■ [0,oo) by 



Proposition 10.4. For every y £ [y~,y^]: the sum W^y- ^y+-^{y) is absolutely convergent. 
Moreover, W^y- ^y+] is continuous wtth respect to pointwise convergence. 

Proof. The compactness of Bui implies that there is a constant _D > so that \dkSj\ < D 
uniformly on Thus, we compute that for y^,y^ G [y~,y~^], 

E E rdkS,ity' + {l-t)y^)dt \yl-yl\<{2r + l)''D J2 \yl-y'^\- 
,„ ,,,, Jo , ,„ 



First of all, this implies that W^[y- y+](j/) is well-defined and converges absolutely for y G 
[y~,y+], because W^[y- .y+j (j/") = 0, by definition, and Efcez-'/H^ jj/fc | < 1 by Theorem 

mm 

Secondly, it is now clear that W^y- y+^ is continuous for pointwise convergence, because 
when y" G [y~,y'^] is a sequence of configurationsconverging pointwise, say to then 
lim„^oo I]jgzd /H \yl '~yT\ = 0' ^ qaite easy to prove, so that lim„_>oo Wyy- y+j (y") = 

The next result is harder to prove: 



Theorem 10.5. We have > 0. Moreover, if a configuration y G [y",?/"*"] is a, 

global minimizer, then W^y- y+^{y) — 0. In particular, W^y- y+^{y^) — .„+] (y^) ~ 0. 

Proof. The proof of this theorem is similar to that of Theorem |4.9| We will sketch it here. 



In fact, we will show that when y^,y^ G [y~,y"'"] and y^ is a global minimizer, then 
^lv~ ,y+](y^^ — "^ly- ,y+](y'^)- Applied to y^ ~ y~ , this shows that W[y~ y+] > 0, whereas 
when applied to y^ ~ y , it shows that W^y- y+]{y) = if y is a global minimizer. 

So let y\y^ G [y-,y+] and suppose that e ~ W^y^ ^y+^{y^) - W^y- ^y+-^{y^) > 0. It 
suffices to show that this implies that y^ is not a global minimizer. To prove this, let 
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(pi, qi), . . . , {pc qc) be principal periods for u and let Bp C be the fundamental domain 
for I'^/Hu, defined by 

Bp:={i£Z'^\Q< {i,pj) < 1 for aU 1 < j < c } . 

Then = Wb,\j^^^ and W„B,{y^) - l^ns,(y') ^ e ■ n\ 

Now define for every n > r + 2 the configuration £ [j/" , by 



j/j^ otherwise. 



Then j/" — is clearly supported in the r-interior of _B„p and it is not too hard to show, 
using the compactness of [y~ ,y'^], the uniform ii-bound on [j/ , y^] and the argument of 
Theorem 14.91 that there is a constant _E > so that 



This means that y^ is not a global minimizer. □ 



One can in fact also prove a variant of Theorem 4.8 that says that if Wyy- y+^{y) = 0, then 
1/ is a global minimizer. Since we do not need this result in this paper, we will not prove it 
here. 

Recall that both [{/~, j/"*"] and [y~ , J/+] are invariant under the forward flow of the negative 
gradient vector field —VW. But a ghost circle is also invariant under the backward flow. 
This implies that, if F is a ghost circle and y G F are the endpoints of a gap in an 

Aubry-Mather set contained in F, then F'*' ' is both forward and backward invariant 
under the negative gradient flow. In order to prove that y~ and y^ axe not the only fixed 
points in y'^], we will now show that W'jy- acts as a Lyapunov function: 

Lemma 10.6. Let y € f'" '^^i and denote hy t '^t the flow of -VW. Then t ^ 
Wyy- y+^ly^ty) is continuously differentiable and 

_ Wiy-,y^]i^ty) = - J2 idM{y)f. 



dt 



Proof. Let us denote, for convenience, ||VW(y)|p :— "Y^ifz^dju {diW{'y))'^. We will begin 
by showing that the function y i-)- ||VM^(j/)||^ is absolutely converger it and continuous on 
[?/~, 3/+]. This is proved by interpolation, as in the proof of Proposition 10.4 More precisely, 
if \diSj'\ < D and \di,kSj \ < C uniformly on [y-,y+], then \diW\ < (2r + ifD, so that 

\\\VW{y')\\^-\\VW{y^)f\ < Yl \d^Wiy') + d,W{y')\-\dMiy')-d.W(y')\ 
< Yl mr + lfD J2 E rd..kS,{ry' + {l-r)y^)dr\-\yl^yl\ 

ieZ-i/H^ \\j~i\\<r\\k-j\\<r ■'^ ' 



Applied to y^ — y~ and y^ = y and combined with Theorem 10.2 this implies that 
||VVy(y)jP is absolutely convergent, because ||VVy(i/~)|p = 0. The continuity for point- 
wise convergence follows from the argument given in Proposition |10.4[ In particular, we 
now know that t i— >■ — ||VIV(^'t?/)|p is continuous, being the composition of two continuous 
functions. 

The next step is to denote y(t) := '^ty, for y £ F'^ and to observe that 

SAyit)) - SM = I' ^S,iy{r))dr = - ^ J' d,S,{y{r)) ■ 9,l^(y(r))dr. 

" ll'=-j||<'" 

Summing this over j we obtain because of the absolute convergence, that 

- = - f \\VW{y{r))fdr. 

Jo 

In particular, because the integrand is continuous, we find that ^ | j_q ^^[j/- (^'tj/) = 
lim,^oi/oi!Viy(y(r))|pdr = -jlVI^(y)lp. " ' □ 
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Proposition 1 10. 4] and Lemma [ 1 . 6| combined now lead to the main result of this section: 
Theorem 10.7. Let [y~ ,y'^] be a gap in the Aubry-Mather set M{x) and let V be a ghost 
circle so that M{x) C F. Then either F'^ ' = F n [y~ ,y'^] consists of global minimizers 
only, or there is at least one stationary point y G F'^ ' that is not a global minimizer. 



10.4 



Proof. By Proposition 
the flow-invariance 

of pfe implies that W^[y- (^ty) = for aU y £ '^^1 and all 



we have that W^[a-,j,+] > 0. When |p[a- ,„+] = 0, then 



t G R. By Proposition [lOe] we then have that ||VlV(y)l|^ = 0. That is, F'^ ' consists of 
stationary points only, and hence by Theorem |10.1| it consists of global minimizers only. 
Because P'^ is compact and y+j is continuous, the other possibility is that 

>i/+] Irlf" f+1 s-ssurnes a positive maximum at some point y £ F'^ '^^^ with y^ <gi y <^ 
y^ . Proposition 10.4 implies that this y is not a global minimizer. It is clearly stationary 
though: if not, then ||VVK(j/)||^ > 0, so that by continuity olt^ 1 1 VM^('I'ty)| we have for 
each t <Q that 

- = - ||VW^(*.y)||^dr > 0. 

Jo 

Because ^ty G P'^ '^^'i this contradicts that y is a maximizer of ,y+] |p[y- y+] • C 

At this moment it is unclear to us whether a gap in an Aubry-Mather set can be foliated by 
stationary points - which therefore all have to be nonrecurrent global minimizers. 



A Twist maps 

Variational monotone recurrence relations do not only arise in statistical mechanics or as 
discretized PDEs: the case of dimension d = 1 is relevant for the theory of twist maps of 
the cylinder. The latter arise for instance in the study of convex billiards and as Poincare 
maps of Hamiltonian systems. In this short descriptive appendix, we will briefly review 
these topics. The informed reader can skip this appendix and we refer to [13] or [9] for more 
detailed proofs of our statements, as well as for a more comprehensive introduction to the 
topic. 

Let us denote by A := E/Z x R the standard cylinder, with coordinates {x mod 1, y) and 
bundle projection tt : A — > R/Z given by {x mod l,y) x mod 1. The lift fr : — >■ R of tt 
to the universal covering spaces sends (x, y) to x. 

Recall that a smooth cylinder map T ; A — >■ A allows for a lift T : R^ — >■ R^, with 
the property that T(x,y) mod (1,0) = T{x mod l,y). This implies that T{x + l,j/) = 
T{x, y) + (n, 0), where n is the degree of T, and moreover that T is unique modulo constants 
of the form (m, 0) with m £ Z. 

Definition A.l. We call a cylinder map T : A — >■ A an exact symplectic positive twist map 
if it satisfies conditions 1, 2 and 3 below. 

1. Degree one: f{x + 1, y) = f{x, y) + (1, 0). 

Condition 1 is true if and only if T is homotopic to the identity map of A. 

2. Exact symplectic: The one-form T*{ydx) — ydx is exact. 

Condition 2 implies that there is a so-called generating function, denoted s : A — >■ R, such 
that T* {ydx) — ydx = ds. We will denote its lift by s : R^ — >■ R. This lift satisfies the identity 
s{x, y) = s{x mod 1, y). In particular, s{x -f- 1, y) = s{x, y). 

Geometrically, condition 2 can be interpreted as follows: one can show that when condi- 
tion 1 holds, then condition 2 is true if and only if T preserves the volume form dy A dx and 
moreover has the property that the volume enclosed by the cycles ^{t) :— {t mod 1, 0) and 
its homotopic image T o 7 is equal to zero. 

Moreover, conditions 1 and 2 hold if and only if T is a so-called Hamiltonian map, i.e. 
T is the time-1 flow of a time-l-periodic Hamiltonian vector field Xu(t) on A. 

To formulate the last condition, let us call T{x,y) = {X{x,y),Y{x,y)). 

3. Positive Twist: The map {x,y) 1-^ (x,X{x,y)) : R^ — ^ R^ is a diffeomorphism. This 
implies that dyX 7^ 0. We require that dyX > 0. 
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Condition 3 says that T twists each fiber {x mod 1} x R C A around the cyhnder A "in 
the positive direction". We will denote the inverse of the map {x,y) i— > {x, X{x,y)) by 
{x,X) ^ {x,y{x,X)). 

In fact, condition 3 allows us to define the function S : — >■ R by S{x, X) :— s{x, y{x, X)). 
The function S is called the generating function of the twist map T. 

The following well-known theorem is crucial in the theory of twist maps. It states that, 
in order to find orbits of exact symplectic positive twist maps of the cylinder, one needs to 
solve a variational monotone recurrence relation in dimension d = 1. For completeness, we 
have included a brief proof of this statement. 

Theorem A. 2. The sequence i {xi,yi) G R^ is an orbit ofT if and only if 

i) For all M, N £ "L with M < N the sequence i t-^ Xi is a stationary point of the finite 
action 

N 

Wm,n{x) ~ ^ Si{x), with Si{x) ■— S{xi-i,Xi), 
for variations of x with fixed endpoints xm and xn ■ 



ii) It holds that yi — —dxS{xi, Xi+i) for all i. 



Moreover, one has that S{x + 1, X + 1) — S{x,X) and dx,xS < 0. 



Proof. Recall the notation T'(x, y) — j/), y(a::, j/)), the diffeomorphism (x, j/) {x,X{x,y)) 

with its inverse {x,X) t-^ {x,y{x, X)) and the definition S{x,X) = s{x,y{x, X)). Then the 
equality T*{ydx) = ydx + ds in the coordinate system {x,y) becomes YdX = ydx + dS in 
the coordinates (x, X), viewing y = y{x, X) and Y = Y{x, y{x, X)) as functions of x and X. 
Writing dS — dxSdx + dxSdX, we thus obtain that 



Y = dxS and y = -d^S 



(A.15) 



Now let i I— >■ {xi,yi) be an arbitrary sequence and define Xi :— X{xi-i,yi-i) and Yi := 
Y{xi-i,yi-i). Then {xi,yi) is an orbit of T if and only if Xi — Xi and yi — Yi for all i. 
According to formula (A.15 1, we have that Yi — Y{xi^i,yi^i) = dxS{xi-i, X{xi-i, i/i-i)) = 
dxS{xi-i,Xi) and yi = —dxS(xi, X{xi,yi)) = —dxS{xi, Xi+i). Thus, Xi = Xi and yi = Yi 
if and only if yt = —dxS{xi, Xi+i) and 



dxS{xi^i,Xi) + dxS{xi,Xi+i) = for aU i G 



(A.16) 



Formula I A. 16 1 holds if and only \ii\-^ Xi is stationary for all the Wn,m defined above. 

Because T has degree one, we have that X{x + 1,1/) = X(x,y) + 1. The function 
y(x,X) is defined implicitly by the relation X(x,y(x, X)) = X, and therefore we see that 
y{x + 1,X + 1) = y{x,X). Thus, the generating function satisfies Six + 1,X + 1) = 
~s{x + 1, y{x + 1, X + 1)) = ~s{x, y{x, X)) = S{x, X). 

Finally, formula (A.15 1 implies that dx,xS — —dxy = —{dyX)~^ < 0. □ 



Examples of twist maps 

Perhaps the most famous example of an exact symplectic twist map is the Chirikov standard 
map. Given a 1-periodic function V — V{x), it is defined as Tv : A — > A by 

Tv{x, y) = {x + y + 2V' (x) mod 1, t/ + 2V'{x)) . 

It turns out that its generating function is S{x,X) := ^{x — X)^ + 2V{x). In other words, 
the variational monotone recurrence relation corresponding to Tv is exactly the Frenkel- 
Kontorova equation in dimension d = 1, given by Si{x) = j{xi-i — Xi)"^ + V{xi). By the 
way, the "standard" is to choose V{x) = cos(27ra;), for some parameter fc > 0. This 
produces the map [x mod 1, y) i— > (a; + j/ — ^ sin(27ra::) mod 1, J/ — ^ sin(27ra;)). 

Another application of the theory of twist maps arises in the context of convex billiards, 
cf. [16| . The configuration space of such a billiard consists of the arclength parameters 
a; G R/Z that describe the position of the billiard ball along the boundary of the billiard 
at the moment of reflection and angles y G (0, vr) measuring the direction of the outgoing 
billiard trajectory with respect to the tangent line to the billiard at x. Then the motion of a 
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Figure 1: An ellipse-shaped convex billiard and part of a billiard trajectory. 



billiard ball is described by an exact symplectic positive twist map T : {xi, yi) i— >■ (sj+i, J/i+i)- 
The variational structure of this problem follows as the rule "angle of incidence = angle of 
reflection" is derived from the variational principle that a billiard ball travels alone "shortest 
paths". The positive twist condition > should be obvious from Figure 1 

Finally, under generic conditions, the Poincare return map of a 2 degree of freedom 
Hamiltonian system near an elliptic equilibrium point is an exact symplectic twist map. In 
this case, the corresponding twist map is actually close to integrable, so that it allows for the 
application of various kinds of perturbation theory. Again, we refer to [ISj for more details. 
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